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ABSTRACT 

We have constructed a sample of 29 close projected quasar pairs where the background quasar 
spectrum reveals absorption from optically thick H I gas associated with the foreground quasar. 
These unique sightlines allow us to study the quasar circumgalactic medium (CGM) in absorp- 
tion and emission simultaneously^ because the background quasar pinpoints large concentrations 
of gas where Lya emission, resulting from quasar-powered fluorescence, resonant Lya scattering, 
and/or cooling radiation, is expected. A sensitive search (la surface-brightness limits of SBLya — 
3 X 10~^^ ergs~^ cm~^ arcsec"^) for diffuse Lya emission in the environments of the foreground (pre- 
dominantly radio-quiet) quasars is conducted using Gemini/GMOS and Keck/LRIS slit spectroscopy. 
We fail to detect large-scale ^ 100 kpc Lya emission, neither at the location of the optically thick 
absorbers nor in the foreground quasar halos, in all cases except a single system. We interpret these 
non-detections as evidence that the gas detected in absorption is shadowed from the quasar UV radi- 
ation due to obscuration effects, which are frequently invoked in unified models of AGN. Small-scale 
Ri. ^ 50 kpc extended Lya nebulosities are detected in 34% of our sample, which are likely the high- 
redshift analogs of the extended emission-line regions (EELRs) commonly observed around low-redshift 
{z < 0.5) quasars. This may be fluorescent recombination radiation from a population of very dense 
clouds with a low covering fraction illuminated by the quasar. We also detect a compact high rest- 
frame equivalent width {Wi^ya > 50A) Lya-emitter with luminosity Li^ya = 2.1 ± 0.32 x lO^^ergs"^ 
at small impact parameter R± = 134 kpc from one foreground quasar, and argue that it is more likely 
to result from quasar-powered fluorescence, than simply be a star-forming galaxy clustered around 
the quasar. Our observations imply that much deeper integrations with upcoming integral- field spec- 
trometers such as MUSE and KCWI will be able to routinely detect a diffuse Lya glow around bright 
quasars on scales R ~ 100 kpc and thus directly image the CGM. 

Subject headings: quasars: absorption lines — galaxies: halos — galaxies: formation — intergalactic 
medium — cosmology: observations 



1. INTRODUCTION 

The ionizing radiation emitted by a luminous quasar 
can, like a flashlight, illuminate hydrogen in its vicin- 
ity, teaching us about the physical properties of gas in 
its surroundings. This is because photoionized hydrogen 
ultimately produces Lya recombination photons, which 
are likely to escape the system. Whether this 'fluores- 
cent' emission arises from the highly ionized Lya forest 
clouds in the intergalactic medium (IGM) which are op- 
tically thin to ionizing photons, from the optically thick 
Lyman Limit Systems (LLSs) and damped Lyman-a sys- 
tems (DLAs) detected as the strongest absorption lines in 
quasar spectra, or from the interstellar or circumgalactic 
medium (CGM) of the quasar host itself, the underlying 
principle is nevertheless the same - a fraction of the en- 
ergy in the quasars' UV continuum is 'focused' into line 
radiation and re-emitted back into space, informing us 
about the physical conditions of the emitting gas. 

The idea of searching for fluorescent Lya emission from 
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the intergalactic medium is not new. iHogan fc Wevm annI 
(1987) first proposed that fluorescent emission might be 
detectable from the optically thin A^hi ^ 17.2, Lya for- 
est clouds, which emit recombination photons because 
they are illu minated by the amb i ent e xtra^alactic UV 
background. ^Gould fc Wein berg (1996) made the cru- 
cial point that emission from the A^hi ^ 17.2 LSSs 
should be much brighter, because these clouds, being 
optically thick to Lyman continuum photons, absorb all 
the ionizing radiation incident upon them. Some ~ 60% 
of these ionizing photons result in Lya recombinations, 
and are 'mirrored' back as Lya photons. Thus at the 
very faintest flux levels, UV background radiation should 
power a ubiquitous fluorescence signal, whereby all of 
the optically th ick g;as in the cosmic web would 'g:low' 
in th e Lya line ("Hogan fc Wevmann"19 8j; iBinette et al.l 
[19931: i Gould fc Weinberg 1996; Cantal upoeTall [2QQ5I : 
"Kollmeier et al.. .2QQ8i ). 

Observationally, a number of increasingly sensitive 
searches for fluorescent radiation powere d by the ambient 
UV background have been perfo r med (iLowenthal et gj] 
IT99QI : iBunker et all IT998I . Il999l : iMarleau et al.l Il999f ). 
the deepest being the ~ 100 hour integration by 
iRauch et all f2008) using the Very Large Telescope 
(VLT). Notwithstanding these efforts, the fluores- 
cence signal has yet to be detected, and given 
independent estimates for the UV background in- 
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tensity (iMeiksin fc W^j ]2QM iBolton et all I2QQ5I : 
iFaucher-Giguere et alJ l2QQ8f ). the expected surface 
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brightness SBLya ^ 1 x 10~"" ergs~" cm~" arcsec 
probably out of reach of current instrumentation. But in 
the regions proximate to a quasar the ionizing flux and re- 
sulting fluorescent surface brightness can be significantly 
enhanced, provided the quasar does not photoevaporate 
the nearby optically thick clouds. Fluorescence in the 
vicinity o f a quasar has been searched for by a number of 
auth ors f Fvnbo et al.' ^199^, 'Francis & Bland-Hawthorn] 
2QQ4l:[FYancis fc McDo n nell 2006; Ad elberger et al. 2006i; 



Cantalupo et al.l 120071 : iRauch et al.l 120081) yielding ei- 



ther null o r questionable det ections. Perhaps most con- 
tentious is lAdelberger et all (p006), who purport to de- 
tect fluorescent emission from a DLA serendipitously de- 
tected in a background quasar spectrum which is 49" 
away from a luminous foreground quasar at z = 2.84. 
Owing to the proximity of the foreground quasar, the 
ionizing flux should be enhanced by a factor of ^ 5700 
over the UVB. However, the large continuum luminosity 
of this source and its implied relatively modest rest-frame 
Lya equivalent width (EW) of 75A suggest that it may 
simply be Lya emission from the galaxy counterpart to 
the DLA, rather t han quasar powe r ed fluo rescence. 

Most recently, I Cantalupo et al.l ()2012l ) performed a 
deep survey for Lya emission around a very bright z = 
2.4 quasar. In their narrow-band imaging, they report 
the detection of ^ 100 sources in their 3.5 Mpc x 3.5 Mpc 
field of view, corresponding to a ~ 5500 comoving Mpc^ 
volume. They further note that a significant subset 
(~ 20%) exhibit Lya EWs far in excess of the max- 
imum predicted value from star-forming regions (i.e. 
Wi^ya > 240 A). These sources are the most convincing 
candidates to date for fluorescing gas in the extended en- 
vironment (several Mpc separation) of, and powered by, 
a quasar. However it still remains unclear whether these 
candidate fluorescent emitters are counterparts to typi- 
cal optically thick absorbers (LLSs and DLAs), and why 
more bright fluorescent emitters were not detected very 
close (< 300 kpc proper) to the quasar. We will address 
both of these questions in this paper. 

Besides illuminating nearby clouds in the IGM, a 
quasar may irradiate gas i n its own host galaxy or 
circumgalactic halo. iReesI ((1988) suggested that cold 
CGM gas in the quasar host illuminated by the quasar 
could be detectable as a n extended "fuzz" of fluorescent 
Lya emission ( s ee als o lAlam fc Miralda-Escudel 120021 : 
iHaiman fc ReesI l2001f ). A number of groups have 
reported the detection of extended Lya emission in the 
vicinity of z ^ 2 — 4 quasars jp iorgovski et al. 1985; 



Hu fc Cow ie 1987; He ckman et al.l ll991b. a: Bremer et al.l 
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19981: Bunker et al." "2002 
Franc is fc Bland-Hawthor^ 
I2OO5I: Francis fc McDonneill 
20061 : Courbin ~al.l 120081 : 
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Hennawi et al.l l200l |North et al. 2012). These efforts, 
some targeted and some serendipitous, reveal a diver- 
sity both in emission level, detection frequency, and 
physical extent, but detailed quantitative interpretation 
is hampered by differing methodologies (deep longslit 
spectroscopy versus narrow-band imaging), sample 
inhomogeneities (radio- loud versus radio-quiet), a broad 
range of redshifts (significant given (1 + z)^ surface 



brightness dimming and strong evolution in the quasar 
population), and ambiguities in detection criteria and 
definitions of the depth of the observations. Neverthe- 
less, the emerging picture is that roughly half of quasars 
from z ^ 2 — 4 exhibit extended emission on scales of 
10 — 50 kpc down to observed frame surface brightness 
levels of SBLya ~ few x 10~^^ergs~^ cm~^ arcsec"^. 
There is suggestive evidence that radio-loud quasars 
have brighter ernission a nd a higher detection frequency 
(jHeckman et al.l 11991b'), although high detection fre- 
quencies have also be en re ported in samples which are 
mostly radio-quiet (IChristensen et al.l 120061 : iNorth e t al.l 
1^12). Finally, this Lya emission may be powered by the 
same mechanism powering the extended emission line 
regions (EELRs) detected aroun d low redshift z < 0.5 
type-I (e.g. Stockton et al. 200 6); iHusemann et al.ll2012h 
and type-II (Greene e t al.l I201lh quasars, traced by 
[O ni] and Balmer lines. 

The brightest cases of Lya nebulosities around 
quasars have an average surface brightness SBLya ~ 
10~^^ ergs~^ cm~^ arcsec"^ {z = 2) over a diameter 
of ~ 100 kpc, corresponding to a total Lya luminos- 
ity of I/Lya ~ 10^^-^ ergs~^. These nebulae bear 
a striking resemblance to the extended Lya emission 
typically obser ved around high -redshift radio galax- 
ies (HzRG; e.g;. iMcCarthv et al.l [19 90: McCarthv 1993' 
van Oiik et al. 1996; Nesvadba et al. 2006; Binette e t al.l 
2006; " Reuland et al. 2007; Villar-Martm et al. 20071 
iMilev fc De Breuckl I2008D as well as the so-called Lya 



'Uoh s' (LABs; e.g. 
Franc is et aH I2OO 



Fvnbo et al.l[T 999l: 'Steidel et al."200C 



iPalunas et al. 2004; Matsuda e t al.l 



2004; Dev et al."'2005"; ^Colbert et al. .2006. : .Nilsson et al.l 
2006; Smith fc Jarvis 2007). The primary difference be- 
tween the nebulae around quasars versus those around 
HzRGs or Lya blobs, is that a strong source of ionizing 
photons (e.g. the quasar) is directly identified, but not 
for the HzRGs or blobs. Obscuration and orientation 
effects, as are often invoked in unifi ed models of AGN 
(see e.g. " Antonuccil [TUl; lElvisI l2000f ) . could be respon- 
sible for this difference. Indeed, the re is ample evidence 
for o bscured AGN in HzRGs (e.g. 'Mil ev fc De Breuckl 
I2OO8I ). and evidence for obscured AGN have been un- 
covered in several of the LABs (|ChaDman et al.l [20041; 
"Basu-Zvch fc S char f "2004"; "Dev et al.'"2005; "Geach et^ 
J007; Barrio et al. 2008; Smith et al. 2008), although 
there are several notable exceptions (Sm ith fc JarvisI 
120071 : [Nilsson et al. 2006: Ouchi et al. 2009). 

Although photoionization by an AGN is perhaps 
the most plausible mechanism powering extended Lya 
nebulae, particularly when the quasar is directly ob- 
served, a litany of other mechanisms have also been put 
forward. Other sources of ionizing photons , such as fast 
radia t ive shocks powered by radio- jets (|McCarthv et al.l 
1 19901 : iHeckm an et al.l Il991bl al) or starburst driven su- 



perwinds ([ Francis et al.l 1200 iTlTaniguchi fc Shioval 120001 : 
iTaniguchi et al.. .200Th have been discussed. It also been 
suggested that X-rays produced via inverse Compton 
scattering of CMB photons (and/or local far- IR pho- 
tons) by relativist ic electrons in radio jets ( Scharf et al.l 
I2OO3I : iSmail et "aTl [2009) . could be the source of ioniza- 
tion. The possibility that extended Lya emission is 
powered by star-formation, either via stellar ionizin g 
photons escaping from galaxies ( Ranch et al.l l2011[ ). 
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spatially extende d star-formation (|Ouchi et al.l I2QQ9I : 
iRauch et al1l2Q12D, or by pure scattering of Lyg photons 
jDiikstra et al." '2006a','c; "Barnes & Haehnelt" "2009", '20T0'; 
[Barnes et al. 2011; Zheng et al. 2011; Steidel et al. 2011; 
l^ijks tra & K ramen i2012), has also been considered. 
In addition to photoionization and scattering, a large 
body of theoretical work has suggested that Lya emis- 
sion nebulae could result from L ya cooling radiation - 
powered by gravitati onal collapse (iHaiman et all 120001: 
' Fardal et al? 2001; Furlanetto et al." "2005"; "Yang et alj 
gQOa [P iikstra et al. 2006b; Diikstra & Loeb 2009f 
Goerdt et al. 2010; D aval et al.l 120101: IBertone fc Schavg 
2010: Faucher-Giguere et al.l 120101 : iFrank et al.l 120121 : 
Rosdahl fc BlaizotI 12012^ " This scenario is particu- 



larly fashionable in light of a plethora of work on 
the so-called 'cold mode' of cosm ological accretion 
(iFardal et an '2001: Katz et al. 2003; Ker es et al 
20091: iDekeffc Birn b oim 2006; Birnboi m et al 



2005 



2007; 



Ocvirk et al'n2008l : iBrooks et all l2009t iDekel et all 



2009*) , whereby a significant fraction of the gas accreting 
onto galaxies from the IGM remains cold T ~ 10^ K, 
and funnels in along large-scale filaments. In the 
absence of significant metal-enrichment, collisionally 
excited Lya is the primary coolant of T ~ 10^ K 
gas; hence the accreting gas could steadily radi- 
ate gravitational potential energy in the Lya line 
as it falls into the halo. While many studies have 
suggested that gravitational co oling radiation could 
powe r the LABs (iFardal et all I2OOII: iFurlanetto et al l 
20051: iDiikstra fc Loebl [20091: IGoerdt et al.l [2010 : 



Faucher-Giguere et al.l 120101 : [Rosdahl fc BlaizQt[ I2OI2 
Frank et al.l |2012[ ). the predictions of these simu- 



lations are uncertain b y orde rs of magnitude (e.g. 
[Furlanetto et al. 2005; Faucher-Giguere et al. 2010; 
Ilosdahl fc Blaizot 2012; Hennawi fc Prochaska 2013) 
because the emissivity of collisionally excited Lya is 
exponentially sensitive to gas temperature. Accurate 
prediction of the temperature requires solving a cou- 
pled radiative transfer and hydrodynamics problem 
which is not curr ently computational feasible (but see 
[Rosdahl fc Blaiio t 2012). 

Thus despite significant observational and theoretical 
efforts, the relationship between Lya nebulae in quasars, 
HzRGs, and LABs is rather confusing, and the physical 
process powering extended Lya emission remains an im- 
portant unsolved problem. Questions also remain about 
Lya emission outside of quasar halos, on the larger scales 
relevant to quasar powered IGM fiuorescence. 

In this fourth paper of the 'Quasars Probing Quasars' 
series, we introduce a novel techniqu43 to tackle the 
important problem of emission from the CGM. We 
search for this emission from close projected quasar 
pairs, which have small angular separation on the sky 
but large line-of-sight separations such that the two 
quasars are physically unassociated. In these unique 
sightlines, absorption in the background ('b/g') quasar 
encodes information about the distribution of circum- 

We note that a similar technique was employed for 
the single serendipitously discovered projected quasar pair by 
(Adelberger et al. 2006), however our approach differs significantly 
in that we use a large sample of quasar pairs to statistically char- 
acterize the properties and distribution of gas in quasar environ- 
ments, which can be used to make predictions for the expected 
emission. 



galactic and intergalactic gas in the vicinity of the fore- 
ground ('f/g') quasar. This technique thus allows us 
to analyze the absorption-line and emission-line prop- 
erties of gas around quasars simultaneously. We per- 
form a systematic, spectroscopic survey for extended 
Lya emission in the vicinity {R± < 600 kpc) of z ^ 
2 quasars, to typical la surface-brightness limits of 
SBLya ^ 3 X 10~^^ ergs~^ cm~^ arcsec"^. The data 
are drawn frorn our s urvey of projected quasar pairs 
([Hennawi et al.|[200 6a1. and we restrict our analysis to 
29 projected pair systems for which the spectrum of the 
b/g quasar shows evidence for an LLS coincident with 
the redshift of the f/g quasar. If the f/g quasar is illu- 
minating this optically thick gas, this is a prime candi- 
date for bright Lya fluorescence and the the b/g quasar 
pinpoints the expected location of the emission. Inde- 
pendently, our observations provide a sensitive search 
for any extended Lya emission in the environment of 
z ^ 2 quasars, such as might be powered by Lya cool- 
ing radiation or any of the other physical mechanism just 
discussed. Our b/g quasar absorption-line measurements 
statistically map out the covering factor of optically thick 
absor ption in the quasar C GM (Hennawi et al. 200^; 
[Hennawi Prochaskal[2007l : iProchaska et al.ll201 2). and 
detailed absorption-line modelin g places constraints on 
the physical properties of this gas ([Prochaska fc Hennawil 
2009). Armed with this knowledge about the properties 
of the CGM, we can make direct predictions for the ex- 
pected fluorescent emission and compare to our observa- 
tions. Similarly, by combining our emission constraints 
from this study with the physical properties of the gas, 
we can directly constrain the heating rate of the CGM 
gas we detect in absorption, which is th e subject of a 
companion paper (Hen nawi fc Prochaskal[2Q^1 3). 

As we will frequently refer to other results from the 
'Quasars Probing Quasars' (QPQ) series, we briefly re- 
view the basic approach and s ummarize the key re- 
sults of each paper. In QPQl ([Hennawi et al.l l2006a[ ) 
we introduced a novel technique to study the physi- 
cal state of gas in the CGM of z o:^ 2 — 3 quasars, 
which are hosted by massive galaxies, using pr ojected 
quasa r pairs (see a lso Crotts 1989: Moller fc Kiaergaardj 
ll992l : lBowen et al.ll2 006: Sc hulze et al.ll2012l : lFarina et al.l 
'2012). Spectroscopic observations of the b/g quasar in 
each pair reveals the nature of the IGM transverse to 
the f/g quasar on scales of a few 10 kpc to several Mpc. 
In QPQl, we searched 149 b/g quasar spectra for opti- 
cally thick absorption in the vicinity of 1.8 < Zfg < 4.0 
luminous f/g quasars, and found a high-incidence of op- 
tically thick gas in the quasar hosts CGM. In QPQ2, we 
compared the statistics of this optically thick absorption 
in b/g sightlines to that observed along the line-of-sight, 
and argued that the f/g quasars emit their ionizing radi- 
ation anisotropcially or intermittently. An echelle spec- 
trum of a projected quasar p air was presented in QPQ3 
([Prochaska fc Hennawil [20091 ) . which resolved the veloc- 
ity field and revealed the physical characteristics (metal- 
licity, dens ity, etc.) of t he absorbing CGM gas. In QPQ5 
(Prochask a et al.ll2012l ). we used a significantly enlarged 
sample of projected quasar pair spectra to characterize 
the CGM of quasars with much improved statistics, fo- 
cusing on the covering factor and EW distributions of 
H I and metal line absorbers. 

This paper is organized as follows. In 0we present a 
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formalism which describes Lya emission powered by pho- 
toionization and scattering. In ^we provide details of 
the quasar pair survey, our spectroscopic observations, 
and describe the selection of the quasar-absorber sam- 
ple that was used to search for extended Lya emission. 
As we will show, a key quantity for interpreting the Lya 
emission properties of the f/g quasars is the covering fac- 
tor of optically thick absorbers measured in background 
sightlines, which we also quantify in this section. A cus- 
tom algorithm was developed to subtract off the quasar 
PSFs and search for faint diffuse extended emission which 
is described in ^ Our results are presented and dis- 
cussed in §5l where we also compare to previous searches 
for Lya fluorescence. We summarize and conclude in 
^ In Appendix A, we quantify the attenuation of Lya 
by dust grains for gas clouds with properties matching 
our absorbers, and conclude that dust extinction is not 
a significant effect. Throughout this paper we use a cos- 
mological model with = 0.24, = 0.76, h = 0.70 
whic h is consistent with t he most recent WMAP5 val- 
ues (|Komatsu et al.|[2009f ). Unless otherwise specified, 
all distances are proper. It is helpful to remember that 
in the chosen cosmology, at redshift z = 2.0, an angular 
separation of = V corresponds to a proper transverse 
separation of R_\_ = 8.7 kpc, and a velocity difference of 
1500 km s~^ corresponds to a radial redshift space dis- 
tance of s = 7.7 Mpc. A typical quasar at 2: = 2.0 with 
an SDSS magnitude of i = 19, has specific luminosity 
at the Lyman edge of log^^o ^j^ll — 30.5 (ergs~^ Hz~^), 
and at a distance of = 100 kpc (11.5^'' projected on 
the sky) the flux of ionizing photons would be a factor 
^uv — 3760 times higher than the ambient extragalac- 
tic UV background. Finally, we use the terms optically 
thick absorbers and LLSs interchangeably, both referring 
to quasar absorption line systems with logA^ni > 17.2, 
making them optically thick at the H I Lyman limit 
(tll > 1). 

2. FORMALISM AND PRELIMINARIES 

This section provides analytical estimates for the sur- 
face brightness of Lya emission from a variety of phys- 
ical processes. We frame the discussion in terms of a 
circumgalactic medium filled with cool clouds, and we 
present results in terms of quantities which can be di- 
rectly measured from absorption-line observations using 
background sightlines. 

2.1. Cool Cloud Model 

Bright emission in the Lya line typically occurs in 
gas at T ~ 10^ K, thus we consider a highly idealized 
model of a population of cool gas clouds in quasar ha- 
los. Several stu dies have previousl y explore d similar 
models (Mo & M iralda-Escudd 119961 : Mailer fc Bullockl 
12004), and others have focused specifically on ex- 
tended Lya emissio n CRees'"1988"; 'Haiman & Rees| |2001|; 
iHaiman et all 120001 : Ipijkstra & Loeb 2009). Nearly ah 
of this work assumes the cool gas is in pressure equi- 
librium with a hot tenuous virialized plasma. Al- 
though there may be some evidence for pressure con- 
finement from our ou r abso rption-line measurements 
(|Prochaska fc Hennawil l2009f ). our conclusions about 
Lya emission do not depend on the existence of a hot 
phase. 



We assume that the cool gas clouds observed in ab- 
sorption are spherical with a single uniform hydrogen 
volume density nn, and that they are spatially uniformly 
distributed throughout a spherical halo of radius R. Al- 
though extremely simple, this model is also relevant, to 
within order unity geometric corrections, to more com- 
plicated distributions which vary with radius, or the fil- 
amentary structures expected in the cold accretion pic- 
ture. In what follows, we often employ averages over the 
projected area of the spherical halo on the sky which we 
denote with angle brackets, e.g. the average line-of-sight 
distance through the halo is (s) = l/Aj sdA = 4i?/3. 
Besides the volume density nn and the size of the halo i?, 
two other parameters completely specify the spatial dis- 
tribution of the gas, namely the hydrogen column density 
of the individual clouds A/'h, and the cloud covering factor 
/c, which is similarly defined to be an average over the 
sphere. We focus on these four parameters {R ,nH ,A^h , 
fc) because they are either directly observable or are 
closely related to absorption-line observations. 

All other quantities of interest can then be expressed in 
terms of these four parameters. For example, of particu- 
lar interest for computing the Lya surface brightness is 
the volume filling factor defined as the ratio of the volume 
occupied by the clouds to the total volume, fy = ricK, 
where ric is the number density of clouds and Vc is their 
volume. The definition of the covering factor is 



(1) 



where ^ 



TicCTc, and <Tc is the cloud cross-sectional 



area. This gives fc = nc(Tc4i?/3 for our case of uniformly 
distributed clouds. Note that in general the covering fac- 
tor can be larger than unity. The corresponding volume 
filling factor can be written as 



h 



and the total mass of cool gas is then 



X 



TTR^fcNa 



X 



(2) 



(3) 



where rup is the mass of the proton and X = 0.76 is the 
hydrogen mass fraction. 

As we will consider emission from objects at high red- 
shift, for which the size of the emitting clouds could be 
unresolved by our instrument, the relevant observable is 
the surface brightness ( erg s~^ cm~^ arcsec"^) averaged 
over the projected area of a sphere as measured from 
Earth. The specific intensity of Lya from our spheri- 
cal halo is given by integrating the equation of radiative 
transfer through the gas distribution 



(4) 



where jLya is the volume emissivity per steradian 
(ergs~^ cm~^ ster~^). In this and the following sections 
we will ignore the attenuation of Lya radiation caused 
by resonant scattering of gas within the halo itself or by 
the IGM, as well as extinction caused by dust. Scat- 
tering within the halo itself would only change the spa- 
tial distribution of the emission, and as we are averaging 
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over apertures on the sky, this should produce smah ef- 
fects. We do however consider the resonant scattering 
of a point source of Lya photons from the quasar CGM 
in §2.31 Scattering by th e IGM is hkely n eghgible for 
the redshifts we consider (jZheng et al.l 1201 1). We esti- 
mate the impact of extinction by dust in the Appendix, 
and demonstrate that it should not result in a significant 
reduction of the Lya emission. 

We can then write the average surface brightness ob- 
served from Earth as 



SB 



1 



Lya 



1 



1 



(1 + Z)4 TtR^ 



jhyadV (5) 



where the factor of (1 + accounts for cosmologi- 
cal surface brightness dimming. The corresponding total 
Lya luminosity is 



(6) 



1.3 X 10^ 



3.0 

SBlyq! 



R 



100 kpc 



10 erg s ^ cm ^ arcsec ^ 



2.2. Lya Fluorescence 



ergs 



For a low density gas, it is possible to separate con- 
tributions to the total Lya intensity into re combina- 
tion and collisionally excited components ( Ch amberlainI 
|1953[ ). Observationally, extended Lya emission could 
have contributions from both mechanisms. Collisional 
excitation of neutral atoms and the subsequent emission 
of Lya can play an important role in the radiative cool- 
ing of T ~ 10^ K gas, and we explore this cooling radi- 
ation from the CGM of quasars in det ail in a compan- 
ion paper (|Hennawi fc Prochaska|[2Q13f ). In what follows 
we consider recombination, which we will also refer to 
as Lya fiuorescence, from a source of ionizing photons. 
A central point source with an dependence is as- 
sumed, but our estimates can easily be generalized to 
other spatial distributions. We consider the two limit- 
ing cases where the gas clouds are either optically thin 
A^Hi < 10^^-^ cm"2 or optically thick Nm > 10^^-^ cm"^ 
to Lyman continuum photons. At intermediate neutral 
columns ^ Nm ^ 10^^*^ cm~^ our simple approxima- 
tions break down and more detailed radiative transfer 
calculations are necessary. In each case, we consider 
the limit where light from the quasar has otherwise been 
unattenuated by resonant scattering or dust. 

2.2.1. Optically Thin 
The recombination emissivity from the cool clouds is 



47r 



(7) 



where Uq is the electron density, is the proton den- 
sity, r^thin = 0-42 is the fraction of recombinations 
which result in a Lya pho t on in the optically thin 
limit (jOsterbrock fc Ferlandl 120061 : IGould fc Weinberg! 
1 19961 ), and the factor of fy accounts for the fact that 
the emitting clouds fill only a fraction of the vol- 
ume. The case A recombination coefficient a a is a 
weak function of temperature T, and we evaluate it 



at T = 10, 000 K where a a = 4.18 x lO-^^cm-^s"^ 
(jOsterbrock fc Ferlandl I2006D . 

The electron and proton densities are determined by 
photoionization equilibrium. 



nmF 



and the photoionization rate F is given by 
1 



(8) 



(9) 



Here ^^ll is the frequency at the Lyman edge, and a^, is 
the hydrogen photoionization cross-section. We assume 
that the quasar's spectral energy distribution obeys the 
power-law form -^ui^i, (z^/z^ll) blue ward of z^ll 

and adopt a slope of ag = 1.57 consistent with the mea- 
surements of EEEilS] (HQQI) . The quasar ionizing lu- 
minosity is then parameterized by Lj^^^^ the luminosity 
at the Lyman ed^e . Throug hout this paper we estimate 
I/j^LL by tying the Te lfer et a l. (2002) power-law spectrum 
to the composite quasar spectrum of Va nden Berk et al] 
(|2001f ), which can then be compared to the SDSS opti- 
cal photometry covering rest-frame UV wavelengths (see 
Appendix A of QPQl for a detailed description of this 
procedure). 

For all cases of interest to us, radiation fields are suf- 
ficiently intense that optically thin gas will always be 
highly-ionized. Under this assumption, the neutral frac- 
tion xm = ^ <C 1, hence ^ nn and ne ~ (1 + Y/ 
2X)nH, where the helium and hydrogen mass fractions 
are Y and X, respectively, and we have assumed all he- 
lium is doubly ionized. 

Putting this all together, the observed Lya surface 
brightness is 

SBL„ = |;^J^a.(l + ^)»„/ciV„ (10) 
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Note that photoionization equilibrium implies that 
jhya FriHi, but siucc the neutral fraction is inversely 
proportional to photoionization rate xri ~ <^a^h(1 + Y/ 
2X)/r, the net effect is that Lya emission is independent 
of the luminosity of the ionizing source in the highly- 
ionized, optically-thin regime, provided that the ionizing 
source is bright enough to keep the gas highly ionized 
and optically thin. This fact can be exploited to obtain 
an expression for the neutral column density averaged 
over the area of the halo, (A^hi), in terms of the observed 
surface brightness and ionizing luminosity: 



(A^Hi) = 2.2 X 10 
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17 
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1030 ergs-iRz-i 



(11) 
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In order for the optically thin approximation to 
be valid, the individual clouds must have A^hi 
10^^"^cm~^. In optically thin photoionization equilib- 
rium, the neutral fraction scales as xm oc so that the 
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largest neutral columns will be attained at the edge 
of the halo. It thus follows that (A/'hi) < max(A/'Hi) and 
hence if (Nm) > 10^^ cm~^ then the clouds are defi- 
nitely optically thick at some locations in the halo and 
the optically thin approximation breaks down. One can 
thus use eqn. (pTj) to compute (A/'hi) and directly test 
whether the optically thin approximation is valid, given 
only the observed Lya surface brightness and the size 
of the nebulae, provided that Lj^^^ can be estimated, as 
in the case of Lya nebulae observed around quasars. If 
(A^Hi) > 10^^-^ cm^ then the clouds cannot be optically 
thin. Note however that the converse is not true — both 
optically thin and optically thick clouds can result in 
(A^Hi) < lO^^-^cm"^ - a smah value of (A^m) thus pro- 
vides no definitive information about whether we are in 
the optically thin or thick regime. 

2.2.2. Optically Thick 

If the clouds are optically thick to ionizing radiation, 
they will no longer emit Lya proportional to the vol- 
ume that they occupy because of self-shielding effects. 
Instead, a thin highly ionized skin will develop around 
each cloud, and nearly all recombinations and resulting 
Lya photons will originate in this skin. The cloud will 
then behave like a special 'mirror', converting a frac- 
tion 77thick = 0.66 of the ionizing photons it intercepts 
into Lya photons emitted at a uniform brightness (i.e 
the specific intensity leav ing the surface is isotropic, see 
iGould fc Weinber"gl Il996f ). Simulations including both 
ionizing radiative transfer and Lya resonant scatter- 
ing (Cantalupo et al. 2005|; iKollmeieT et al. 2008) have 
shown, that for the case of a quasar illuminating an op- 
tically thick cloud, the true surface brightness is reduced 
from the 'mirror' value by a factor of 2 — 4. This re- 
duction occurs because illumination and radiative trans- 
fer effects redistribute photons over a wide solid angle. 
We introduce a geometric reduction factor /gm to ac- 
count for this reduction. For the case of a transversely 
illuminated spherical 'mirror', the resulting 'half- moon' 
illumination pattern res ults in a /gm = 4/(37r) = 0.42 
(|Adelberger et al.ll2006[ ). We set %^ = 0.5, which is 
consistent wi th the range seen in radiative transfer sim- 
ulations (Ca ntalupo et al.l 120051 : Ikollmeier et a l. 2008). 
At a distance R from a quasar, a mirror will emit a sur- 
face brightness 
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where $ (phots ^ cm ^) is the ionizing photon number 
fiux 



(13) 



The total observed flux from the mirror is simply pro- 
portional to the area of the mirror times this SBLya- 

Now consider the case of a population of cool optically 
thick clouds, hence also mirrors, but which are spatially 



unresolved by our instrument. In this regime we can 
write the Lya volume emissivity as 



^thick/i^Lya ^ 



(14) 



Note that there is no geometric reduction factor in this 
case. Conservation of photons requires that a fraction 
^thick of the ionizing photons emerge as Lya photons, 
and the emissivity must be isotropic. Again integrating 
the equation of radiative transfer and averaging over area 
(eqns. |4]and[5]), we obtain 
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In the expression above it is understood that fc is capped 
at /c = 1, corresponding to the case where all ionizing 
photons are absorbed. 

Thus in the optically thick limit, the emission depends 
only weakly on the amount of cool gas through the cover- 
ing factor of the clouds, and scales linearly with the ioniz- 
ing photon flux. This contrasts strongly with the highly 
ionized optically thin limit, where the surface brightness 
depends on the amount of cold gas through the combi- 
nation /c^hA/^h, but not on the intensity of the ionizing 
source. 

2.3. Lya Emission from Scattering 

In this subsection we compute the surface brightness of 
extended Lya emission, produced via resonant scattering 
of a central source of Lya photons (i.e. the quasar) by the 
neutral gas in the CGM. By analogy with the optically 
thick fluorescence, the scattering surface brightness can 
be written as 



(16) 



where fc now represents the covering factor of gas which 
is optically thick in the Lya transition (i.e. A/'hi ^ 
lO^^cm"^). This will in general be much higher than 
the covering factor of clouds optically thick at the Ly- 
man limit considered in ^2.2.21 The quantity <^>Lya is the 
flux of ionizing photons emitted close enough to the Lya 
resonance to be scattered by gas in motion around the 
quasar 



L 



Lya 



(17) 



where we have introduced the rest-frame absorption 
equivalent width Wi^ya , and implicitly assumed that it is 
a very weak function of radius, as we demonstrate next. 

The volume averaged neutral column density through 
our gas distribution is (A^m) = /cA^h(^hi), where {xm) 
is the volume averaged neutral fraction determined from 
photoionization equilibrium (eqn. [8]). These relations im- 
ply an average Lya optical depth 
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(18) 



where we have introduced the 6-value describing turbu- 
lent and thermal motions of the absorbing gas. On the 
flat portion of the curve of growth 1 < TLv g 5. IQf i th e 
corresponding equivalent width is (see e.g. iDrainell2QlL) 



Wi^ya - 0.99 



50 km s ^ , 



Vln((TLyc.)/ln2) 
2.1 



A, 

(19) 

SO that equivalent widths of Wi^ya — lA are expected. 
We can finally write 



SBLyc. = 1.2 X 10 



WLy 
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(20) 



where it is again understood that /c < 1 in eqn. (|2Q]) . 

Because equivalent width is such a weak function of 
optical depth on the saturated part of the curve of 
growth (eqn. [T9|) . WLya is approximately independent 
of the ionizing luminosity of the quasar as well as the 
details of the gas distribution. Thus the scattering sur- 
face brightness scales as SBLya cx /c^Lya, completely 
analogous to optically thick recombinations for which 
SBLya (X fc^ (eqn. [T5]) . However, scattering will typ- 
ically be about ~ 100 times fainter than optically thick 
recombinations if the covering factors are the same be- 
cause 77thick*^*/*^*LyQ! ^ 100, which rcsults from the as- 
sumption of Wi^ya — 1.0 A, and from the relative amount 
of ionizing photons to Lya photons emitted by a typi- 
cal quasar. One expects this scattering emission to be 
important, therefore, only in regimes where the gas is 
optically thin at the Lyman limit. Indeed, the surface 
brightness of scattering (eqn. [2Q|) and optically thin fiuo- 
rescence (eqn. [TQ|) can be comparable, but have very dif- 
ferent dependencies on gas distribution and quasar lumi- 
nosity. Throughout this work we estimate Lj^^^^ by tying 
the c omposite quasar spectrum of ( Vand en Berk et al. 
I200lf ) to the SDSS optical photometry (covering rest- 
frame UV wavelengths) and computing the peak of the 
Lya emission line. 

2.4. An Illustrative Example 

There are two important points that can be made 
about the surface brightness of extended Lya emission, 
and its implications for the supply of cool gas. First, op- 
tically thin and optically thick fluorescence can result in 
comparable levels of emission, as can pure resonant scat- 
tering of a central source of Lya photons. To complicate 
matters further, cooling radiation, that is collisional ex- 
citation of Lya, can also lead to comparable emission. In 
general, the only way to determine which mechanism is at 
work is to use diagnostics from additional emission lines. 
For instance, the detection of strong extended Balmer 
line emission would be provide compelling evidence for 
recombinations, and effectively rule out scattering and 
cooling. But given that Balmer lines are redshifted into 
the near-IR for z > 2, detecting even Ha would be a 



daunting observational task because of the low expected 
surface brightness. Furthermore, even if one is certain 
that recombination is the physical mechanism, it is still 
typically impossible to ascertain the ionization state of 
the gas, and hence whether it is optically thin or thic^El. 

The second and related point is that measurements 
of the Lya surface brightness alone cannot be used to 
determine the cool gas mass. For recombinations, the 
surface brightness scales as SBLya cx /c^h^h if the gas 
is optically thin (eqn. [10]), or as SBLya cx $/c if it is 
optically thick (eqn. [15]). Whereas, for scattering the 
scaling is SBLya cx ^Lyafc (eqn. [20]). The total cool gas 
mass scales with the different combination Mc oc fc^u 
(eqn. [3]). Even if we have perfect knowledge of the mech- 
anism (recombinations versus scattering) and the ion- 
ization state (optically thin or optically thick), the cool 
gas mass is not uniquely specified by the Lya emission 
alone. Note that the situation is even worse for Lya 
cooling radiation because of an e xponential dependence 
on an unknown gas temperature (jHennawi fc Prochaskal 
[2011 . Thus even in the best-case scenario where Balmer 
line observations point to recombinations, one will not 
know whether the gas is optically thin or thick, and hence 
which dependence on physical parameters to use. 

The following worked example illustrates the afore- 
mentioned issues. Consider a relatively bright z = 2 
quasar with apparent magni tude of i = 17.0 . Us ing the 
procedure in Appendix A of iHennawi et al.l (|200 6a). we 
determine the specific luminosity to be, log^^g ^z^ll ~ 



31.3 at the Lyman limit, and log^^g ^ 



10 -^Z^Lya 



32.0 in 



Lya. Now consider three distinct cool gas distribu- 
tions, one optically thin to ionizing photons with param- 
eters {R,nu,Nu/fc) = (100 kpc, 0.15 cm-^, lO^^cm-^, 
1.0), another optically thick to ionizing photons with 
{R.nu^NuJc) = (100 kpc, 100cm-^ lO^^cm-^, 0.01), 
and a third optically thin to ionizing photons but op- 
tically thick to Lya photons with (i^, nn, A^h, /c) = 
(100 kpc, 30 cm-^, 10^^ cm-^, 2.0). These very dif- 
ferent cool gas distributions result in nearly identi- 
cal average Lya surface brightness SBLya ^ 1.1 — 
1.2 X 10~^^ ergs~^ cm~^ arcsec"^ over an = 100 kpc 
halo. For the first two models the emission is powered by 
photoionization, whereas for the third it is pure Lya scat- 
tering. But note also that for the first distribution, Lya 
scattering contributes at very nearly the same surface 
brightness level ~ 1.1 x 10~^^ ergs~^ cm~^ arcsec"^ as 
optically thin photoionization, implying that disentan- 
gling the two mechanisms via say detection of Balmer 
lines would be challenging. These sizes and surface 
brightness levels are in line with the brightest Lya blobs 
surveyed by Matsuda et alj (|2004 1201 ID (rescaled to 
2: = 2) as well as the extended Lya nebulae detected 
around z - 2 HzRGs (e.g. Vihar-Martm et alj I2007D 
and radio- loud quasars (e.g. iHeckman et al.til991bl ). For 
the first optically thin model the total cool gas mass 
is Mc = 3.3 X 10^^ M0, for the second optically thick 
model it is Mc = 2.5 x 10^ M©, and for the Lya scat- 
tering case it is Mc = 6.6 x 10^ M©. Thus variations 
of ~ 5000 in the total gas mass can all lead to compa- 

^ For extreme ratios of SBLya to Ljy-^-^ , it is possible to rule out 
optically thin clouds via eqn. However, this requires that the 

ionizing flux can be measured, which is possible for quasars but 
not HzRGs or LABs. 
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rable levels of extended Lya emission, and even if one 
were sure that photoionization was powering the emis- 
sion (from Balmer lines) and/or detected the source of 
ionizing photons (an unobscured quasar), there is no 
way to determine whether the gas is optically thin or 
thick. Plugging the observed observed surface bright- 
ness SBLya, distance and ionizing luminosity L^^^l 
into eqn. (pTj) for the area averaged neutral column gives 
(^Hi) = 1-3 X lO^^cm"^, and thus we cannot distinguish 
an optically thin from an optically thick scenario. 

2.5. Why Quasar Pairs? 

In the previous subsection we saw that a variety of 
physical processes and a wide range of gas distributions 
can all lead to comparable Lya emission. Hence, the 
great value of quasar pairs is that absorption-line obser- 
vations of b/g quasars can provide crucial information 
about the physical properties of the gas in the CGM 
of the f/g quasar. Indeed, all the physical parame- 
ters (i^ , nn , A^H , /c) can be determined from absorp- 
tion line observations, provided that the covering factor 
is high enough for the clouds to be occasionally inter- 
cepted by b/g sight lines. For example, we map out the 
covering factor of optically thick absorption as a func- 
tion of impact parameter (see Figure [2|) in ^3.41 (see 
also iHennawi et all l2QQ6at iHennawi" ProchasyMrzl : 
iProchaska et al.l l2Q12V and find a high covering factor 
of cool, op tically t hick ^ as around all bright quasars at 
2. In IProchaska fc Hennawil (|2QQ9 ). we conducted 
detailed absorption line modeling of a single projected 
pair and determined that the column densities of the 
absorbing clouds are typically Nn — 10^^ cm~^ and vol- 
ume densities nu 0.1 — lcm~^, although it is unclear 
whether this single system is representative. Armed with 
this knowledge of the gas distribution, one can then inter- 
pret the physical implications of the presence or absence 
of extended line emission. 

Note however that the properties of gas with a very 
low covering factor fc ~ 0.01 are much more challeng- 
ing to measure, since thousands of background sources 
would be required to map it out. Thus a crucial caveat 
in comparing absorption line observations to emission 
line observations is that the gas responsible for the de- 
tected emission may arise from clouds with a much lower 
covering factor than the gas dominating the absorption. 
Nevertheless, one can still make quantitative statements 
about the expected emission from the component of gas 
which is detected in absorption, which will have the high- 
est covering factor. This is the approach adopted in this 
work. 

But even if quasar pairs are useful for measuring the 
properties of the quasar CGM, why use quasar pairs, as 
opposed to individual isolated quasars, to searching for 
diffuse emission lines? Consider the following four sce- 
narios which are illustrated schematically in Figure [TJ 
a) The absorbers/emitters are large resolved clouds at a 
'random' location in the quasar CGM (FigurefTJi). In this 
case, the background sightlines pinpoint the locations of 
gas and thereby determine where to l ook for emission. 
This a pproach is akin to that taken by lAdelberger et al.l 
(|2006l) , who purported to detect quasar fluorescence from 
an optically thick absorber detected in a b/g sight line. 
For this case of large resolved absorbers, the relevant SB 
estimate is the mirror approximation (see eqn. [12]). b) 



Another possibility is that the clouds are small and hence 
unresolved by our instrument, but that the gas resides in 
large resolved filamentary complexes which cover a frac- 
tion of the quasar halo (Figure [TJd). For this case, the b/g 
sightline intersects a filament by construction, and iden- 
tifies the location where emission is expected via longslit 
observations of the quasar pair. Note that the average 
SBs computed in ^ would underestimate the true SB in 
this case, since they include a reduction by the covering 
factor /c, due to regions devoid of gas. c) It could also 
be that the clouds are small and hence unresolved, but 
some quasars have cool gas complexes and others do not. 
This bimodality in gas supply could result in the 'hits' 
and 'misses' which give rise to the observed covering fac- 
tor (Figure dh). In this case, the quasar pairs are helpful 
since the background sightline identifies which objects 
to search for emission around, and the average SBs that 
we computed in ^are again underestimated because the 
covering factor includes objects with no gas. d) Finally, 
the clouds could be small and hence unresolved, and dis- 
tributed rather uniformly in the quasar halo (Figure [l]i). 
In this case the background sightline does not provide us 
with useful information about where to search for emis- 
sion. The average SB appropriately accounts for the di- 
lution by the covering factor /c, and this is the SB that 
we expect to detect at any spatial position along the slit. 
For this case, a search for diffuse emission around iso- 
lated individual quasars, i.e. not residing in pairs, would 
be just as effective. 

3. QUASAR PAIR OBSERVATIONS 

Our goal is to conduct a sensitive search for ex- 
tended Lya emission from the CGM and IGM around 
f/g quasars in quasar pairs, and to relate these emission 
constraints to our knowledge of the distribution of gas 
from absorption-line measurements of b/g sightlines. We 
concentrate on a unique sample of projected quasar pairs 
useful for exploring the CGM/IGM surrounding high- 2; 
quasars and, by inference, the massive galaxies and dark 
matter halos which host them. We search for extended 
Lya emission around 29 z ~ 2 quasars, which are drawn 
from a parent sample of 68 projected quasar pairs. Be- 
fore delving into details, we briefly summarize the criteria 
used to select these 29 quasar pairs: 

• the CGM/IGM around the f/g quasar exhibits 
strong and optically thick H I absorption in the 
b/g quasar sightline. 

• the transverse separation of the pair and 
magnitude of the f/g quasar imply a flu- 
orescent mirror surface brightness SBLya > 
5 X 10""*^^ ergs""*^ cm~^ arcsec"^ (see eqn. [T2]). 

• a deep spectroscopic observation of the f/g quasar 
was obtained covering Lya , allowing us to perform 
a sensitive search for Lya emission on scales of tens 
to hundreds of kpc. 

• the two quasars are separated by velocity > 
2000 km s~^, to ensure that they are indeed pro- 
jected and not physically associated. 

Although our search for extended Lya emission is fo- 
cused on sightlines showing strong absorbers, we char- 
acterize the covering factor of optically thick absorbers 
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Figure 1. Four scenarios for the spatial distribution of CGM gas in quasar halos which illustrate the usefulness of projected quasar pair 
observations, a) The absorbers/emitters are large resolved clouds at random locations in the quasar CGM. The background sightlines 
pinpoint the locations of optically thick clouds and thereby deter mine where to look for emission. For such large spatially resolved clouds, 
the relevant SB estimate is the mirror approximation (see eqn. I12|) . b) The clouds are small and hence spatially unresolved by our 
instrument, but gas resides in large resolved filamentary complexes which cover a fraction of the quasar halo. For longslit observations of 
quasar pairs, the b/g sightline intersects the filament by construction, again pinpointing where emission is expected. The average SBs that 
we computed in jSl underestimate the true SB of the filament, since they include a reduction by the covering factor /c, due to regions 
devoid of gas. c) The clouds are small and unresolved, but some quasars have cool gas complexes and others do not, and this bimodality in 
the gas supply is responsible for the variations that result in the observed covering factor. Quasar pairs are helpful because the background 
sightline identifies objects with large gas supply and hence where to search for emission. The average SBs that we computed in ^are again 
underestimated because the covering factor includes objects with no gas. d) The clouds are small, unresolved, and distributed uniformly 
in the quasar halo. The background sightline does not provide useful information about where to search for emission. The average SB 
estimates in ^2] ctppropriately account for the dilution by the covering factor /c, and this is the surface brightness that we expect to detect 
at any spatial position along the slit. For this case, a search for diffuse emission around isolated individual quasars, i.e. not residing in 
pairs, would be just as effective. 



using the full unbiased parent sample of 68 projected 
quasar pairs. This analysis demonstrates that the CGM 
surrounding z ^ 2 quasars has a high > 50% covering 
fract ion of optically thick gas to radius R ^ 200 kpc (see 
also iHennawi et al] l2006ar iHennawT Prochaskdl2007l : 
iProchaska et al.ll2012[ ). 

There are several reasons why we restrict our search 
for emission to only pair sightlines showing strong ab- 
sorbers. First, in the context of fluorescence from the 
IGM, powered either by the UVB or boosted by a nearby 
quasar, emission from optically thick clouds will always 
dominate over optically thin gas because the density of 
gas in the optically thin Ly a forest is just too low (e.g. 
iGould fc Weinberd Uool see also eqn. [TQ|) . Although 



this may not hold true in the CGM, where gas densi- 
ties can be higher, note that the the mirror SB from 
a spatially resolved optically thick cloud is independent 
of the gas properties (eqn. [12]), and the same holds for 
emission from a population of optically thick clouds pro- 
vided the covering factor has been measured (eqn. [T5|) . 
In contrast, optically thin emission depends on the com- 
bination /c^H^H which is typically unknown. The up- 
shot is that a non-detection of emission from gas which 
is optically thin could result either from /c^h^h be- 
ing too low, or because the gas is not illuminated by 
the quasar. Whereas, interpretation of non-detection of 
emission from optically thick gas is unambiguous, and 



10 



HENNAWI & PROCHASKA 



implies that the gas is not ihuminated by the quasai0. 

Note however, that given the high covering factor of op- 
ticahy thick gas in the quasar CGM, statistically, there 
is a high probability that any long-slit observation of a 
projected quasar pair will intersect optically thick mate- 
rial, even if the b/g sightline does not show absorption. 
We refer the reader to the discussion in ^2.51 and Fig- 
ure [U which illustrates, nevertheless, why having a b/g 
sightline showing optically thick absorption can be ad- 
vantageous. 

In the following subsections we describe the details of 
the quasar pair survey, the spectroscopic observations, 
the definition of the quasar-absorber sample, and present 
our measurement of the covering factor of optically thick 
absorbers from the parent sample. 

3.1. Quasar Pair Survey 

Modern spectroscopic surveys select against close pairs 
of quasars because of fiber collisions. For the Sloan Digi- 
tal Sky Survey (SDSS), the finite size of optical fibers pre- 
cludes discovery of pairs with separation < 55^' Q Thus, 
to find pairs with sub-arcminute separations, additional 
follow-up spectroscopy is required both to spectroscopi- 
cally confirm quasar pair candidates, and to obtain spec- 
tra of sufficient quality to search for absorption line sys- 
tems. 

In an ongoing survey, close quasar pair candidates are 
selected fro m a photometric quasar catalog (B ow et al.. 
[20TTL [20121 ). and are confirmed via spectroscopy on 4m 
class telescopes including: the 3.5m telescope at Apache 
Point Observatory (APO), the Mayall 4m telescope at 
Kitt Peak National Observatory (KPNO), the Multi- 
ple Mirror 6.5m Telescope, and the Calar Alto Obser- 
vatory (CAHA) 3.5m telescope. Our continuing effor t 
to discover qua sar pairs is de scribed in l Hen nawi' f200l), 
iHennawi et all jlOOGb), and iHennawi et al. (2010). To 
date about 350 pairs of quasars have been uncovered with 
impact parameter R±^ < 300 kpc and z^g > 1.^. Pro- 
jected pair sight lines were then observed with Keck and 
Gemini to obtain 'science quality' high signal-to-noise 
ratio (S/N) moderate resolution spectra. A subset of 
higher resolution spectra at echellette and echelle reso- 
lution have also been obtained, but these are not used 
in this work. This observational program has several 
science goals: measure the s mall-scal e quasar cluster- 
ing!; of quasars (|Hennawi et al. 2006b; Mvers et al.l 120081 : 
IHennawi et al.ll2Q10l : iShen et al...2009.,) . to measure small 
scale transverse Lya forest correlations (Rorai et al. in 
prep), to characterize the transverse proximity effect 
(Hennawi et al., in prep.), and to use the b/g sightline to 
characterize the c ircumg al actic medium of the f/^ quasar 
(IHenna wi et~all l2006at IHennawi Prochaskal boOTi : 
IProchask afc Hennawil 120091 : iProchaska et al.l I2012D . 
which is relevant to the goals of this paper. 

^ Another possibility is extinction of Lyo; by dust, but we show 
in Appendix A that this effect is not significant. 

^ An exception to this rule exists for a fraction (~ 30%) of 
the area of the SDSS spectroscopic survey covered by overlapping 
plates. Because the same area of sky was observed spectroscopi- 
cally on more than one occasion, the effects of fiber collisions are 
reduced. 

^ The lower limit on redshift is motivated by the ability to 
detect redshifted Lyo; absorption above the atmospheric cutoff 
A > 3200 A. 



3.2. Spectroscopic Observations 

High S/N ratio, moderate resolution slit spectra of 
~ 100 quasars were obtained with Keck and Gemini 
in observing runs spanning from 2004 until 2008. All 
quasars observed have z > 1.6, which is the lower limit 
for detecting Lya set by the atmospheric cutoff. About 
two-thirds of the targeted pairs consist of projected pairs 
of quasars {Av > 2500km s~^) at different redshifts; the 
rest were physically associated binary quasars. 

For the Keck observations, we u sed the Low Re solu- 
tion Imaging Spectrograph (LRIS; lOke et al]ll995f ). ei- 
ther in longslit mode targeting only the two quasars in 
the pair, or in multi-slit mode with custom designed slit- 
masks, which allowed placement of slits on other known 
quasars or quasar candidates in the field. Typically the 
slitmask or slit was rotated such that both quasars in 
the close pair could be observed on the same slit. LRIS 
is a double spectrograph with two arms giving simulta- 
neous coverage of the near-UV and red. We used the 
D460 dichroic with the 1200 lines mm~^ grism blazed 
at 3400 A on the blue side, resulting in wavelength cov- 
erage of ~ 3300 — 4200 A. The dispersion of this grism 
IS 0.50 A per pixel and our 1" slits give a resolution of 
FWHMc:^ 150km s~^. These data provide the coverage 
of Lya at Zfg for all of our pairs. On the red side we typ- 
ically used the R600/7500 or R600/10000 gratings with 
a tilt chosen to cover the Mg II emission line at the f/g 
quasar redshift, useful for determining accurate systemic 
redshifts of the quasars. Occasionally the R1200/5000 
grating was also used to give additional bluer wavelength 
coverage. The higher dispersion, better sensitivity, and 
extended coverage in the red provided high signal-to- 
noise ratio spectra of the Mg II emission line and also 
enabled a more sensitive search for metal-line absorption 
in the b/g quasar (see Prochaska et al. in prep.). Some 
of our older data also used the lower-resolution 300/5000 
grating on the red-side covering the wavelength range 
4700 - 10, 000 A. About half of our LRIS observations 
were taken after the atmospheric dispersion corrector was 
installed, which reduced slit-losses (for point sources) in 
the UV. 

The Gemini data were take n with the Gemin i Multi- 
Object Spectrograph (GMOS; "Hook et al."'2004) on the 
Gemini North facility. We used the B1200_G5301 grating 
which has 1200 lines mm~^ and is blazed at 5300 A. The 
detector was binned in the spectral direction resulting in 
a pixel size of 0.47 A, and the f slit corresponds to a 
FWHM 0^ 125km s"\ The slit was rotated so that both 
quasars in a pair could be observed simultaneously. The 
wavelength center depended on the redshift of the quasar 
pair being observed. But most of the the targets consid- 
ered here are at z ~ 2.2 — 2.5, so the grating was typically 
centered at 4500 A, giving coverage from 3750 — 5225 A. 
The Gemini CCD has two gaps in the spectral direction, 
corresponding to 9 A at our resolution. The wavelength 
center was thus dithered by 15-50 A between exposures to 
obtain full wavelength coverage in the gaps. The Gem- 
ini observations were conducted over three classical runs 
during UT 2004 April 21-23, UV 2004 November 16-18, 
and UT 2005 March 13-16 (Program IDs: GN-2004A-C- 
5, GN-2004B-C-4, GN-2005A-C-9, GN-2005A-DD-4). 

Total exposure times ranged from 600 — 21800s, for 
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the Keck and Gemini observations, depending on the 
magnitudes of the targets. This was typicahy broken 
up into several individual exposures with exposure times 
of 600- 1800s depending on the total planned exposure 
time. The motivation for the shorter (total) integrations 
was to conduct a survey to build up a statistical sample 
of b/g absorption line spectra at small impact parameter 
from the f/g quasars. The goal of the longer integrations 
was to obtain high-quality spectra of a subset of quasar 
pairs and to search for low SB Lya emission near the 
f/g quasar. The S/N ratio in Lya forest region of the 
b/g quasar varies considerably, but it is almost always 
S/N> 5 per pixel for the data we consider here. 

3.3. Defining the Quasar-Absorber Sample 

In QPQ5 and Prochaska eFall (HoH) (QPQ6) we an- 
alyze the statistical properties of the CGM near quasars 
using our full dataset of 'science quality' projected- 
pair spectra, wh ich combines spectr a from th e SDSS 
([Abazajian et al. 2005) and BOSS (jAhn et al.l 12012ft 
spectroscopic surveys with higher quality data from 8m- 
class telescopes. Here we consider a subset of this full 
dataset constituting 68 sightlines, and survey the spec- 
tra for optically thick absorption coincident with the f/g 
quasar redshift. We first describe the additional criteria 
applied to the full dataset to arrive at the 68 sightline 
pa rent sample, and then summarize the procedure used 
in iProchaska et al.l ()2012f ) to identify optically thick ab- 
sorbers coincident with the f/g quasars. 

The present study concerns a sensitive search for ex- 
tended Lya emission, which may coincide with the loca- 
tion of an optically thick absorber detected in a neigh- 
boring background sightline. To this end, we restrict our 
parent sample to only the subset of quasars which have 
an implied SBLya > 5 x 10~^^ ergs~^ cm~^ arcsec"^ 
given by the mirror approximation in eqn. ([T2]) . where 
the radius R is set to be the impact parameter R± of the 
projected quasar pair. This surface brightness limit is 
chosen because it approximately matches the median 2a 
surface brightness limit that our longslit spectroscopic 
observations achieve (see the next section for detailed 
discussion of surface brightness hmits). 

The only other criteria imposed are that we must have 
obtained a spectroscopic observation of the quasar pair 
with Gemini or Keck covering Lya at the f/g quasar 
redshift, and we require that the velocity difference be- 
tween the f/g and b/g quasar be > 2000 km s~^ to ensure 
that the two quasar are not physically associated. This 
minimum velocity difference is many times larger than 
both the expected line-of-sight velocity dispersion in a 
z ^ 2 quasar environments ~ 300 km s~^, and the typi- 
cal redshift errors ~ 270 - 660km s"^ (see below). We 
further culled the sample by excluding pairs whose b/g 
quasar exhibits strong broad absorption lines (BAL), as 
evidenced by large C IV or Lya equivalent widths (EWs). 
Mild BALs were also excluded if the BAL absorption 
clearly coincided with the velocity window about Lya at 
the f/g quasar redshift. 

The remaining quasar pairs were visually inspected and 
searched for significant Lya absorption within a velocity 
window of \Av\ = 1500km s~^ about the f/g quasar red- 
shift. This velocity window is chosen to accommodate 
physically associated absorbers with extreme kinemat- 



ics 1000km s~^, as well as to bracket uncertainties in 
the f/g quasar systemic redshifts. Quasar redshifts de- 
termined from the rest-frame ultraviolet emission lines 
(redshifted into the optical at z ~ 2) can differ by up to 
one thousand kilometers per second from the systemic 
frame, because of outflowing /inflowing material in the 
broad line regions of quasa rs (Gaskeh 1982; Tvtler & Fan] 
19921: iV anden Berk et all [2001 : Richards et al. 2002b^ 
Shen et al. 200^. Systemic redshifts are estimated by 
combining the line-centering procedure used in QPQl 
with the recipe in Shen et al. (2007) for combining mea- 
surements from different emission lines. The resulting 
typical 1(7 redshift uncertainties using this technique are 
in the range 270 — 770km s~^ depending on which 
emission lines are used. This analysis considered all avail- 
able spectra on the f/g quasar including the public SDSS 
and BOSS datasets. We refer the reader to QPQ6 for 
further details. 

The Lya transition saturates at A^hi — 10^^ cm~^, and 
between column densities of Nm 10^^~^^cm~^ the 
curve of growth is flat, and Wi^ya is a poor proxy for hy- 
drogen column density. In principle, additional leverage 
may be provided by simultaneously fitting higher order 
Lyman transitions, or the detection of Lyman limit ab- 
sorption, which requires A^hi > lO^^'^cm"^. However, 
nearly all of our sightlines are at redshifts z < 2.6, for 
which the Lyman limit is below the atmospheric cutoffl, 
and so constructing a complete and fully certifiable sam- 
ple of optically thick absorbers with A'hi > 10^^-^ cm~^ 
is extremely challenging given our wavelength coverage. 
These challenges are exacerbated by noise in the spectra 
and the line blending which inevitably occurs at moder- 
ate resolution (FWHM 150km s"^). 

We identified the strongest absorption line in the 
± 1500km s~^ velocity window about the f/g quasar red- 
shift, and measured the Lya equivalent width Wi^ya- All 
systems with Wi^ya > 1 A were flagged for Voigt profile 
fitting of the Lya transition, to estimate the H I col- 
umn density. For reference, a Lya equivalent threshold 
of Wi^ya > 1.5 A {Wi^ya > 2 A) corresponds to column 
densities of roughly A'hi > 10^^ cm-^ (A^m > 10^^ cm-^) 
if the absorption is dominated by a single component; but 
this need not be the case. 

The H I analysis was complemented by a search for 
metal-lines at the f/g quasar redshift in the clean con- 
tinuum region redward of the Lya forest of the b/g 
quasar. The narrow metal-lines provide a more accu- 
rate redshift for the absorption line system (than ob- 
tainable from Lya) and, if present, they can help dis- 
tinguish optically thick absorbers from blended Lya for- 
est lines. We focused on the stronges t low-ion transi- 
tions commonly observed in DLAs (e.g. ^ Prochaska et al.l 
[200Th : Si II A1260, 1304, 1526, O I A1302, C II A1334, 
Al II A1670, Fe II A1608, 2382, 2600, Mg II A2796,2803; 
and the strong high- ionization transitions commonly seen 
in LLSs: C IV A1548, 1550 and Si IV A1393, 1402. 
Any systems with secure metal-line absorption were also 
flagged for Voigt profile fitting. Again, this analysis oc- 
casionally benefited from the public datasets of SDSS 
and BOSS whose broad wavelength coverage comple- 
ments the more limited coverage of our large-aperture 

^ In fact, we rarely even cover Ly/3. 



Table 1 

Properties of Projected Quasar Pairs Sample Surveyed for Diffuse Lya Emission 



Quasar Pair 


B/G Quasar 




F/G Quasar 




^bg 










2^ 

_a_s 




Radio^ 


9^ 




R\ 




SBjnir'^ 


SBi^k 


4x 
exp 


Inst.™ 






logio^^Si 


SDSS J0203-0002 


J020341.34-000235, 


.1 


J020340.40-000201. 


.8 


1, 


.984 


1, 


.735 


272 


1, 


.7335 


-219 


< 0.98 


36. 


.1 


316 


30.4 


1.31 


0.75 


1200 


LRIS 


2.6 + 0. 


.2 


18.9 ± 0, 


.3 


SDSS J0239-0106 


J023946.45-010644, 


.2 


J023946.43-010640. 


,5 


3 


.129 


2 


.292 


272 


2 


.3024 


988 


< 0.98 


3, 


.7 


32 


30.1 


36.14 


0.14 


7200 


LRIS 


10.3 + 0. 


,3 


20.3 ± 0, 


,2 


SDSS J0250— 0047 


J025039. 82— 004749, 


.6 


J025038. 68— 004739. 


,2 


2 


.447 


1, 


.851 


662 


1, 


.8547 


429 


< 0.93 


20. 


,0 


175 


31.0 


14.97 


0.21 


2400 


LRIS 


5.9 ± 0. 


,2 


19.7 ± 0, 


,4 


SDSS J0750+2724 


J075008.27+272404, 


.5 


J075009.25+272405. 


,2 


1. 


.802 


1, 


.771 


272 


1, 


.7691 


-248 


< 0.99 


13. 


,1 


115 


30.0 


4.46 


0.28 


2400 


LRIS 


2.3 ± 0. 


,2 


< 18.9 




SDSS J0752+4011 


J075259. 14+401 118, 


,2 


J075259.81+401128. 


.2 


2, 


.121 


1, 


.883 


272 


1, 


.8830 


-21 


< 0.98 


12. 


.6 


110 


29.9 


2.87 


0.28 


1200 


LRIS 


3.8 ± 0. 


.2 


18.8 ± 0, 


.3 


SDSS J0800+3542 


J080048. 74+354231, 


.3 


J080049. 89+354249. 


.6 


2, 


.066 


1, 


.981 


272 


1, 


.9824 


179 


< 0.97 


23. 


.1 


201 


30.6 


3.82 


0.31 


3139 


LRIS 


2.0 ± 0. 


.1 


18.7 ± 0, 


,4 


SDSS J0814+3250 


J081419. 58+325018, 


.7 


J081420. 37+325016. 


.1 


2, 


.213 


2, 


.173 


792 


2, 


.1790 


574 


< 0.92 


10. 


,3 


88 


30.3 


9.06 


0.32 


12600 


GMOS 


2.3 ± 0. 


,2 


18.8 ± 0, 


,3 


SDSS J0837+3837 


J083757. 13+383722, 


.4 


J083757. 91+383727. 


.1 


2, 


.251 


2, 


.059 


272 


2, 


.0612 


199 


8.61 


10. 


,3 


89 


30.5 


14.60 


0.39 


4530 


LRIS 


2.5 ± 0. 


,1 


18.6 ± 0, 


.3 


SDSS J0841+3921 


J084159. 26+392140, 


.0 


J084158. 47+392121. 


,0 


2, 


.214 


2, 


.041 


272 


2, 


.0475 


654 


< 1.02 


21. 


,1 


183 


30.4 


2.59 


0.23 


1800 


LRIS 


3.9 ± 0. 


,2 


19.2 ± 0, 


.2 


SDSS J0856+1158 


J085656. 05+115802, 


.7 


J085655. 75+115802. 


.0 


2, 


.483 


1, 


.767 


272 


1, 


.7641 


-299 


10.89 


4, 


.5 


40 


30.6 


153.06 


0.15 


4800 


LRIS 


16.0 ± 0. 


.2 


21.0 ± 0, 


.2 


SDSS J0924+3929 


J092416. 72+392914, 


.6 


J092417.65+392920. 


,3 


2 


.080 


1, 


.885 


272 


1, 


.8868 


170 


< 0.97 


12. 


,2 


106 


30.3 


7.42 


0.37 


1800 


LRIS 


1.5 ± 0. 


,1 


< 18.8 




SDSS J0938+5317 


J093804. 84+531743, 


.1 


J093804. 22+531743. 


.9 


2 


.323 


2, 


.063 


272 


2, 


.0675 


412 


< 0.93 


5, 


.6 


49 


29.9 


12.72 


0.29 


1800 


LRIS 


1.6 ± 0. 


.1 


< 18.8 




SDSS J1006+4804 


J100627. 10+480429, 


,9 


J100627. 47+480420. 


,0 


2 


.591 


2 


.297 


272 


2 


.3015 


367 


< 0.94 


10. 


,6 


90 


30.0 


3.48 


0.23 


1800 


LRIS 


4.1 ± 0, 


,2 


19.1 + 0, 


,2 


bDbb J l(Jz5+58z(J 


J 102554.77+582017, 


,0 


J 102553.47+582012. 


.0 


2, 


.567 


1, 


.956 


662 


1, 


.9534 


-287 


< 0.88 


11. 


,4 


99 


30.0 


4.25 


0.37 


1500 


LRIS 


2.4 ± 0, 


,2 


19.1 + 0, 


.3 


SDSS J1041+5630 


J104129.27+563023, 


.5 


J104121. 90+563001. 


.3 


2, 


.264 


2, 


.054 


272 


2, 


.0589 


468 


< 1.03 


65. 


.0 


562 


31.0 


1.08 


0.38 


900 


LRIS 


1.9 + 0, 


,1 


< 18.8 




SDSS J1045+4351 


J104506.39+435115, 


.3 


J104508.88+435118. 


.2 


3, 


.011 


2, 


.428 


272 


2, 


.4366 


726 


< 0.89 


27. 


.1 


229 


30.4 


1.08 


0.36 


1460 


LRIS 


1.7 + 0, 


,1 


< 18.8 




SDSS J1145+0322 


J114546.22+032251, 


.9 


J114546. 54+032236. 


,7 


2, 


.011 


1, 


.765 


272 


1, 


.7644 


-86 


< 1.02 


15. 


,9 


139 


29.9 


2.40 


0.62 


800 


LRIS 


1.5 + 0, 


,2 


< 18.9 




SDSS J1204+0221^ 


J120416.69+022111, 


.0 


J120417.47+022104, 


,7 


2, 


.528 


2, 


.436 


44 


2, 


.4403 


375 


< 0.94 


13. 


,3 


112 


30.2 


2.72 


0.18 


7321 


GMOS 


4.9 + 0, 


,2 


19.6 + 0, 


.2 


SDSS J1306+6158 


J130603.55+615835, 


.2 


J130605. 19+615823. 


.7 


2, 


.175 


2, 


.109 


272 


2, 


.1074 


-136 




16. 


.3 


141 


30.2 


2.82 


0.24 


3600 


LRIS 


10.7 + 0. 


.2 


20.2 + 0, 


.2 


SDSS J1358+2738 


J135849.54+273756, 


.9 


J135849.71+273806. 


,9 


2 


.127 


1, 


.900 


272 


1, 


.8930 


-684 


< 0.92 


10. 


,2 


89 


29.9 


4.65 


0.31 


2400 


LRIS 


2.1 + 0. 


,1 


< 18.8 




SDSS J1413+2715 


J141337.96+271511, 


.0 


J141337.18+271517. 


.1 


1, 


.965 


1, 


.689 


272 


1, 


.6885 


-94 


< 0.91 


12. 


.0 


105 


29.6 


2.06 


0.94 


1860 


LRIS 


1.7 + 0. 


.1 


< 18.8 




SDSS J1420+1603 


J142054.92+160342, 


.9 


J142054.42+160333. 


.3 


2 


.057 


2, 


.020 


272 


2 


.0224 


273 


< 0.96 


12. 


.0 


104 


30.6 


13.99 


0.58 


1050 


LRIS 


5.1+0. 


,2 


< 19.5 




SDSS J1427-0121 


J142758. 74-012136, 


.2 


J142758. 89-012130. 


,4 


2, 


.353 


2, 


.281 


272 


2, 


.2821 


100 


30.63 


6, 


.2 


53 


30.6 


42.30 


0.14 


21800 


GMOS 


2.8 + 0. 


,1 


< 18.8 




SDSS J1431+5727 


J143109. 22+572726, 


.4 


J143109.67+572729. 


.0 


2, 


.063 


1, 


.679 


272 


1, 


.6835 


499 


< 0.95 


4, 


.5 


39 


29.6 


17.64 


0.35 


2400 


LRIS 


2.1 + 0. 


.2 


19.0 + 0, 


.3 


SDSS J1442+0137 


J144231. 91+013734, 


.8 


J144232. 92+013730. 


.4 


2, 


.272 


1, 


.807 


272 


1, 


.8068 


8 


< 0.97 


15. 


.7 


137 


30.3 


5.87 


0.30 


1800 


LRIS 


3.2 + 0. 


.2 


19.1 + 0, 


.3 


SDSS J1444+3113 


J144427.96+311314, 


.0 


J144429.40+311321. 


.2 


1, 


.795 


1, 


.734 


272 


1, 


.7370 


294 


< 0.79 


19. 


.8 


173 


31.2 


28.64 


1.60 


600 


LRIS 


1.5 + 0. 


.1 


< 18.8 




SDSSJ1508+3635 


J150812.80+363530, 


.3 


J150814.06+363529. 


,4 


2 


.105 


1, 


.848 


272 


1, 


.8514 


347 


< 1.00 


15. 


,2 


133 


30.7 


14.45 


0.34 


1800 


LRIS 


2.0 + 0. 


,1 


< 18.8 




SDSS J1553+1921 


J155325.89+192137, 


.7 


J155325. 61+192141. 


,0 


2 


.098 


2, 


.016 


662 


2 


.0130 


-326 


< 0.97 


5, 


.1 


44 


29.7 


9.92 


0.27 


1605 


LRIS 


7.7 + 0. 


,2 


20.1 + 0, 


.2 


SDSSJ1605+5113 


J160546.67+511322, 


.9 


J160547.61+511330. 


,5 


1, 


.844 


1, 


.785 


662 


1, 


.7817 


-342 


< 0.98 


11. 


,6 


102 


30.4 


12.24 


1.08 


600 


LRIS 


1.8 + 0. 


2 


< 18.8 





Note. — 

Redshifts of b/g quasar. 
Redshifts of f/g quasar. 
Error on f/g quasar redshift (km s"-*^). 

Absorber redshift. Determined from metal- line transitions when present; otherwise, it corresponds to the approximate centroid of Lya. 
Velocity difference between absorber and f/g quasar (km s~^). 

Peak 20cm radio flux Fpeak (mJy) of f/g quasars which had a match in the source catalog from the VLA/FIRST survey (Becker et al. 1995). For f/g quasars without FIRST matches we list the 5cr upper 
mit based on the sky-rms coverage maps. One source labeled — was not covered by FIRST. 
Angular separation between f/g and b/g quasar (arcsec). 
Impact parameter at f/g quasar redshift. 
Specific luminosity at the Lyman limit (ergs~^ Hz~^). 

Predicted mirror surface brightness (10~^^ erg s~^ cm~^ arcsec"^) according to eqn. (12). 
Observed la surface brightness limit (10~-^^ erg s~^ cm~^ arcsec"^). 
Exposure time (seconds). 
Rest-frame Lyaequivalent width (A). 
Neutral hydrogen column density (cm~^). 

For SDSSJ 1204++0221 we use the f/g quasar redshift and error from near-IR spectroscopy published in Prochaska Hennawi (2009). 
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Figure 2. Properties of projected quasar pair sample. Middle: 
Distribution of impact parameter and f/g quasar redshift for the 
68 projected pairs in our our parent sample. Filled blue circles 
are 29 sightlines with an optically thick absorber within |A^'| = 
1500km s-i of the f/g quasar redshift, which we search for Lyo; 
emission (see Table 1). The black open circles are sightlines which 
are ambiguous, and the green open circles are sightlines classified 
as optically thin. Bottom: Blue squares show the covering factor of 
optically thick absorbers in logarithmically spaced bins. Vertical 
error bars are Poisson. Because some of the ambiguous points 
(black open circles in middle panel) may also be optically thick, 
our covering factor is a conservative lower-limit, as denoted by 
arrows on the blue points. The lower red error bars indicate the 
expected covering factor within a | Ai;| = 1500km s~^ interval from 
the background abundance of LLSs ( Ribaudo et al.| [2011'). with 
the range given by the current precision of these measurements. 
The covering factor of optically thick absorbers around quasars is 
high > 50% for scales R±^ < 200 kpc, and represents a significant 
enhancement over the cosmologic al ba ckground. Top: Distribution 
of predicted mirror SBmir (eqn. (|12|) for the 29 quasars searched 
for emission. Note the general RJ^ trend. The dashed line is the 

median value of the la SB limits, SBicr of our sample^ T he red 

arrow indicates the effective la SB limit of the Ranch e t al.l (|2008l ) 
92 hour integration, where their values are rescaled from 2; ~ 3.1 to 
the median redshift of our sample z = 2 (i.e. (l-\-z)^ dimming) and 
also rescaled to match our 700 km s~^ velocity aper ture. The green 
arrow indicates the average SB threshold given the IMatsuda et all 
|2004) definition of a LAB (fi^ya > 0.7 x lO^^^erg s"^ cm"^ over 
a 16 arcsec^ isophote) also rescaled to z = 2. 



spectroscopy. 

Given the equivalent width for Lya absorption, our 
Voigt profile fits for the A/'hi, and the presence/absence 
of low-ion metal absorption, objects were classified into 
three categories: optically thick, ambiguous, or optically 
thin. Objects which show obvious damping wings or 
strong (EW > 0.3 A) low- ion metal absorption are clas- 
sified as optically thick. For cases for which metal lines 
are weak, not covered by our spectral coverage, or sig- 
nificantly blended with the Lya forest of the b/g quasar, 
an object is classified as optically thick only if it has 
Wi^ya^ 1.7A. Table 1 lists relevant quantities for the 
final sample of 29 objects classified as optically thick, 
around which we search for extended Lya emission. 

The false positives and completeness of our optically 
thick sample are sources of concern. Line-blending, in 
particular, can significantly depress the continuum near 
the Lya profile biasing the I^Lya and the estimated col- 
umn density high. Based on the analysis of our current 
sample, and comparisons of objects observed at moderate 
and echellette/echelle resolution (see e.g. Hennawi et al.l 
|2006a), we estimate that false-positives, that is non 
optically-thick objects, in our sample is lower than 10%. 
This low contamination follows from our relatively con- 
servative criteria for defining a system as optically thick. 
Specifically, there are 3 cases for which we are less con- 
fident in the optically thick designations (see Table 1; 
J1041+5630, J1045+4351, J1444+3113) but we never- 
theless include them in the following analysis. Removing 
them would make no difference to our conclusions aside 
from reducing the sample size by 10%. Note that sat- 
urated Lya forest lines with A^hi < 10^^-^ cm~^ are es- 
sentially indistinguishable from optically thick absorbers 
with low-metallicity. As such, we expect the false posi- 
tive rate to be smaller than our incompleteness, and we 
thus consider our estimate of the covering factor of op- 
tically thick gas (see Figure [2j) to be conservative lower 
limits. 

Finally, given that there is suggestive evidence that 
radio-loud quasars have brighter Lya emission neb- 
ulae and a higher d etection frequ ency than radio- 
quiet quasars (Hec kman et al.|[l991bl ), and that HzRGs 
at 2: ~ 2 typically exhibit bright (SBLya ~ 
IQ-iQ erg s~^ cm~^ arcs ec~^) la rge-scal e 100 kpc) 

Lya nebulae (e.^. McCarthv et a l. 1990; Mc Carthvl[T99l 
van Oiik et al.lll996l: iNesvadba e t al. 2006; Binet te et all 



20061: i Reulan d et al.l 120071 : iVillar-Martm et al.l 120071 : 



Milev fc De Br eucH 120081 ) . we consider the radio prop' 



erties of our f/g quasars. The Faint Images of the Radio 
Sky at Twenty cm survey (FIRST; Becker et al. 1995|) 
used the Very Large Array (VLA) to produce a map 
of the 20 cm (1.4 GHz) sky with a beam size of 5/^4 
and an rms sensitivity of about 0.15 mJy beam~^. The 
survey covers the same 10,000 deg^ sky region covered 
by the SDSS imaging, has a typical detection threshold 
of 1 mJ y, and an a strometric accuracy of 0.05". Fol- 
lowing (jivezic et al.l [20 02), we match our f/g quasars 
to the FIRST catalog with a matching radius of 1.5". 
From our parent sample of 68 projected pairs, only 59 
sources were covered by the FIRST imaging footprint, 
and of these, 4/59 have matches in the FIRST cata- 
log (Fpeak ^ liiiJy)^ which is consi stent with roughl y 
~ 10% of quasars being radio- loud (|Ivezic et al.ll2002[ ). 
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This is to be expected since our f/g quasars are typi- 
cal SDSS quasars. Visual inspection of the FIRST im- 
ages indicates that none show evidence of a complex 
radio morphology (i.e. core- jet , core-lobe, and double - 
lobe; Magliocchetti et al. "1998"; "McMahon et al.' '2002). 
Among the 29 f/g quasars with coincident optically thick 
absorption, 28 sources were covered by the FIRST foot- 
print, and only 3/28 have matches in the FIRST catalog, 
again consistent with a ~ 10% radio loud fraction. The 
radio fluxes for the three sources detected in the FIRST 
survey are listed in in Table 1 ; for the sources which were 
not detected we list 5cr upper limits based on the FIRST 
sky-rms coverage maps. Given that only ~ 10% of our 
f/g quasars have radio counterparts > 1 mJy, both in 
our parent sample and in the sample with coincident ab- 
sorbers, we conclude that quasars which we study here 
are predominantly radio-quiet. 

3.4. The Covering Factor of Optically Thick Absorbers 

The distribution of f/g quasar redshifts and transverse 
separations for all 68 projected quasar pair sightlines 
in our parent sample are shown in the scatter plot in 
Figure O The 29 black filled symbols indicate sight- 
lines which have an optically thick absorber within an 
interval 5v = [-1500, 1500] km s"^ of the f/g quasar 
redshift (the median absolute offset of these systems is 
\5v\ = 256kms~ ), which are objects that we search 
around for extended Lya emission. Open black circles in- 
dicate sight fines classified as ambiguous, and green sight- 
lines which are optically thin. Note that the closest pairs 
are predominantly at redshift z ~ 2, because this is where 
quasar selec tion, and hence quasar pair selection, is most 
efficient (seelRichards et al.ll2002at iHennawi et al .112006131 : 
iBovv et al.ll201lL I2012D . 

It is clear from Figure [2] that the covering factor of 
absorbers is very high for small separations: at least 28 
out of 48 sightlines with R± < 200 kpc have signatures 
of optically thick gas coincident with the f/g quasar, im- 
plying a > 58% covering fraction . This high c overing 
factor was previously noted in Hennawi et al. (2006a) 
based on a much smaller pair sample (11 sightlines with 
i^perp < 200 kpc) which are a subset of our current par- 
ent sample (see also QPQ5). Note again that because 
of the significant incompleteness discussed in §3.3[ this 
high covering factor should be interpreted as a conserva- 
tive lower limit, albeit subject to Poisson uncertainty. 

The distribution of points in Figure [21 or equivalently 
the dependence of covering factor on impact parameter, 
can be used to measure the quasar-absorber correlation 
function, which is equivalently the spatial profile of the 
absorbing clouds around the f/g quasar. In QPQ2, we an- 
alyzed the much smaller sample of pairs from QPQl and 
determined a correlation length of ro = 9.2t.\[^ Mpc 
(comoving) assuming a power law correlation function, 
^ (X r~^, with 7 = 1.6. The high covering factor of ab- 
sorbers quantified by the correlation function, implies a 
comparably high incidence of proximate absorbers along 
the hne-of-sight; however, this absorption is not observed 
in our f/g quasars, or in other isolated quasars. This 
anisotropic clustering of absorbers around quasars in- 
dicates that the transverse direction is less likely to be 
illuminated by ionizing photons than the line-of-sight, ei- 
ther because of anisotropic emission which could result 



if the f/g quasar is obscured from some vantage points, 
or because the quasar emits radiation intermittently, on 
a timescale comparable to the transverse light-crossing 
time. We will return to this important inference in the 
context of the search for Lya fiuorescence in ^ 

4. DATA REDUCTION AND PSF SUBTRACTION 

Our goal is to conduct a sensitive search for extended 
Lya emission in slit spectra in the CGM region surround- 
ing the f/g quasar in quasar pairs. Our approach is to 
subtract off the point spread function (PSF) of the the 
quasars and other objects (galaxies or stars), which may 
have serendipitously landed on the slit, and search for 
resolved Lya emission which is inconsistent with these 
spectral PSFs. In the next sections, we describe the de- 
tails of our data reduction and spectral PSF subtraction 
algorithm, our procedure for ffux calibration and deter- 
mining surface brightness limits, and our simulation us- 
ing synthetic sources to determine the surface brightness 
levels we can visually detect. 

4.1. Data Reduction and PSF Subtraction Algorithm 

The Gemini/GMOS observations were conducted in 
longslit mode, whereas the LRIS observations include a 
mix of longslit and multislit exposures. All data were 
reduced using the LowRedux pipelin43, which is a pub- 
licly available collection of custom codes written in the 
Interactive Data Language (IDL) for reducing slit spec- 
troscopy. Individual exposures are processed using stan- 
dard techniques, namely they are overscan and bias sub- 
tracted and ffat fielded. Flat fielding is performed in 
two steps, using both a pixel fiat to correct for pixel-to- 
pixel sensitivity variations, as well as spectroscopic illu- 
mination fiats (taken with the sky in twilight) to fiatten 
the larger scale illumination pattern arising from either 
instrument optics or imperfections in the slits. Cosmic 
rays and bad pixels are identified and masked in multiple 
steps. First, a sharpness detection algorithm is run on 
each image to identify and mask features smaller than 
the spectral/spatial PSF. Then further downstream in 
the reduction procedure, outliers from our models of the 
sky and the object are masked (see below). Wavelength 
solutions are determined from low order polynomial fits 
to arc lamp spectra, and then a wavelength map is ob- 
tained by tracing the spatial trajectory of arc lines across 
each slit. This wavelength map allows us to model the 
sky counts as a function of wavelength and slit position 
witho ut needing to rectify the original data (e.g. iKelsonl 
[2003h . 

Our method for spectral PSF subtraction employs a 
novel custom algorithm. We treat sky and PSF subtrac- 
tion as a coupled problem to obtain a multi-component 
model of the image counts. This model is composed of 
a sum of 2-d 'basis functions', which consists of a sky 
background and a PSF model for each object on the slit. 
This is typically just the two quasars (f/g and b/g), but 
in some cases it includes additional objects which landed 
on the slit in the region of interest. By construction, 
the sky-background has a fiat spatial profile because our 
slits are fiattened by the slit illumination function. For 
the object models, we first identify objects in an initial 
sky-subtracted image, and trace their trajectory across 

http:/ /www. ucolick.org/~xavier/LowRedux 
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the detector. We then extract a 1-d spectrum, normahze 
these sky-subtracted images by the total extracted flux, 
and fit a B-spline profile to the normalized spatial light 
profile of each object relative to the position of its trace. 
This non-parametric object profile has the flexibility to 
vary as a function of wavelength, and corrections to the 
initial trace are simultaneously determined and applied. 
Given this set of 2-d basis functions, i.e. the flat sky and 
the object model profiles, we then minimize chi-squared 
for the best set of spectral B-spline coefficients (i.e. break 
points are spaced in wavelength) which are the spectral 
amplitudes of each basis component of the 2-d model. In 
other words 



iVpix 

E 



(DATAi - MODELi 



(21) 



where the sum is taken over all A^pix pixels in the image, 
'DATA' is the image, 'MODEL' is a linear combination 
of 2-d basis functions multiplied by B-spline spectral am- 
plitudes, and cr is a model of the noise in the spectrum, 
i.e. cr^ = SKY + OBJECTS + READNOISE. The result 
of this procedure are then full 2-d models of the sky- 
background (SKY), ah object spectra (OBJECTS), and 
the noise (cr^). We then use this model SKY to update 
the sky-subtraction, the individual object profiles are re- 
fit and the basis functions updated, and chi-square fitting 
is repeated. We iterate this procedure of object profile 
fitting and subsequent chi-squared modeling four times 
until we arrive at our final models for the sky background, 
the 2-d spectrum of each object, and the noise. 

Each exposure of a given target is modeled according 
to the above procedure which then allows us to subtract 
the sky and the object profiles from each individual im- 
age. These images are registered to a common frame by 
applying integer pixel shifts (to avoid correlating errors), 
and are then combined to form final 2-d stacked sky- 
subtracted and sky-and-PSF-subtracted images. The in- 
dividual 2-d frames are optimally weighted by the (S/N)^ 
of their extracted 1-d spectra (using the b/g quasar), 
and bad pixels identified via sharpness filtering, as out- 
liers in the chi-square fitting, or by sigma-clipping of the 
image stack, are masked. The final result of our data 
analysis are three images: 1) an optimally weighted av- 
erage of (DATA - SKY), henceforth the 'stacked sky- 
subtracted image', 2) an optimally weighted average of 
(DATA-SKY-OBJECTS), henceforth the 'stacked sky- 
and-PSF-subtracted image', and 3) the noise model for 
these images <j^. The final noise map is propagated from 
the individual noise model images taking into account 
weighting and pixel masking entirely self-consistently. 

Extended Lya emission will be manifest as residual 
flux in our 2-d sky-and-PSF-subtracted images which is 
inconsistent with being noise. We thus search for Lya 
emission by defining a x image in analogy with eqn. (|2T]) 
above, but using the stacked images and correspond- 
ing propagated noise instead of the individual frames. 
This procedure allows us to visually assess the statisti- 
cal significance of any putative emission feature. Fig- 
ure [5] presents these images for the 29 objects in our 
optically thick sample (see Table 1). The middle pan- 
els show Xsky+PSF = (DATA - SKY - OBJECTS)/cr for 
the stacked sky-and-PSF-subtracted images. Recall that 



an analogous has been minimized by our modeling of 
each image. Thus if our model is an accurate descrip- 
tion of the data, the distribution of pixel values in the 
Xsky+PSF should be a Gaussian with unit variance. The 
lower panel of each image shows Xsky = (DATA — SKY) / 
cr; the numerator is the the stacked sky-subtracted (but 
not PSF-subtracted) image. The upper panels show 
smoothed maps which are helpful for identifying ex- 
tended emission. Specifically, the smoothed images are 
given by 

_ CONVOL[DATA - SKY - OBJECTS] 

CONVOL[a] ' ^ ' 

where the CONVOL operation denotes smoothing of 
the stacked images with a symmetric Gaussian kernel 
(same spatial and spectral widths) with FWHMsmth = 
235 km s ^ (dispersion (Jsmth = 100 km s"^). For LRIS 
this FWHMsmth corresponds to 5.2-6.1 pixels, or 1.4-1.5 
times the spectral resolution element, and 1.4 — 1.6'^ spa- 
tiall}El. For CMOS a FWHMsmth = 235 km s"^ cor- 
responds to 6.6-7.0 pixels or 1.4-1.5 times the spectral 
resolution, or 1.9 — 2.0^' spatially. 

Note that smoothing correlates pixels in the smoothed 
image, and furthermore the ratio of smoothed image to 
smoothed noise will no longer obey Gaussian statistics. 
Nevertheless, Xsmth proves to be a useful tool for identi- 
fying extended emission. 

4.2. Flux Calibration and Surface Brightness Limits 

The following describes our procedure to flux cali- 
brate our spectra and determine surface brightness lim- 
its. Standard star spectra were not typically taken im- 
mediately before/ after our quasar pair observations. In- 
stead, we construct a model sensitivity function of the 
Keck LRIS-B spectrograph by first fitting an aggregate 
of standard star spectra taken at different slit positions, 
which span the range of wavelength coverage available 
with a multi-slit setup. A similar procedure was em- 
ployed for Gemini CMOS using standard star spectra 
taken at a variety of central wavelengths. We apply 
this archived sensitivity function to the b/g quasar spec- 
trum, and then integrate the flux- calibrated 1-d spec- 
trum against the SDSS ii-band filter curve. The sensi- 
tivity function is then rescaled to yield the correct SDSS 
ii-band photometry. A comparison of our 1-d spectra 
flux calibrated in this way to the SDSS spectra typically 
show agreement to within ~ 10 — 20% in the region of 
wavelength overlap. Given that SDSS spectra have spec - 
trophotometric errors of - 10% (|Abazai lan "d~Zll2005h . 
that ix-band photometric errors (used for our renormal- 
ization) are about 7% for a typical source with :^ 21, 
and that quasar variability over 5 year timescales 
(between when the SDSS photometry and our spectra 
were taken) could r esult in a ~ 10% fluctuation (e.g. 
iSchmidt et aII[20Toh . we consider ~ 20% agreement be- 
tween our spectra and the SDSS spectra to be reasonable. 

This procedure of rescaling the sensitivity functions 
to match photometry is effective for point source flux- 
calibration, however a subtlety arises in relation to cali- 
brating extended emission. The point source counts will 

The range of values results from the different f/g quasar red- 
shifts, and hence different observed frame wavelength of Lya. 
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be reduced because of improper object centering due to 
guiding errors or inaccuracies in slit /mask alignment, 
and because a fraction of the spatial PSF lies outside 
the slit area. These point-source slit-losses do not, how- 
ever, reduce the amount of spatially resolved emission, 
which is of interest here. Thus, renormalizing to point- 
source photometry will tend to over-estimate our surface 
brightness detection limit and underestimate our sensi- 
tivity. Hence, our procedure is to apply the rescaled sen- 
sitivity functions (based on point source photometry) to 
our 2-d images, but reduce them by a geometric slit-loss 
factor so that we properly treat extended emission. 

To compute the slit-losses we use the measured spatial 
FWHM to determine the fraction of light going through 
our 1.0'^ slits, but we do not model centering errors. A 
useful check on our flux calibration procedure is to exam- 
ine any residual variation in our renormalized sensitivity 
functions after the effects of the atmospheric extinction 
(i.e. due to airmass) and slit-losses are removed. These 
variations in the spectroscopic zero-point could be due 
either to transparency variations or to other systematic 
flux calibration errors. We find a 35% relative variation 
(Icr) about the mean zero-point, which is at a level that 
can be plausibly explained by transparency variations in 
clear but non-photometric conditions. If these variations 
in zero-point are indeed due to transparency fluctuations, 
then our procedure of renormalizing to the photometry 
effectively takes them out. But because we have no quan- 
titative information about the transparency during our 
observations, we conservatively assume a relative error 
of ~ 35% on our surface brightness limits. 

Application of the wavelength dependent sensitivity 
function converts our stacked sky-and-PSF-subtracted 
images from units of electrons into erg s~^ cm" -2 A-i. We 
then tile these calibrated images with a uniform grid of 
aperture centers within ±3000 km s~^ of the f/g quasar's 
redshifted Lya line. Regions within ±7.0^' of the f/g 
quasar trace on either side are excluded to avoid poten- 
tial contamination from Lya emission. Regions within 
±3000 km s"^ of the b/g quasar's Lya and also within 
±7.0^' of the b/g quasar trace are also excluded, to sim- 
ilarly avoid Lya contamination (this is relevant if the 
redshift difference of the pair Sv < 3000 km s"^). The 
spacing between these aperture centers is 1.0'^ spatially 
and Av = 400 km s~^ spectrally. We then compute 
SB limits in windows of 700kms~^ x 1.0^' about each 
aperture center, which corresponds to an aperture of 
700kms~^x l.Oarcsec^ on the sky because we always 
used a 1.0^' slit. The 700 km s~^ spectral width is mo- 
tivated by the observed kinematics of the metal-line ab- 
sorbers associ ated with the optically thick gas we detect 
in absorption (jProchaska fc Hennawil[2Q09L Prochaska et 
al. in prep.). If this gas emits Lya, we expect the emis- 
sion kinematics to have a comparable velocity width (in 
the absence of significant resonant broadening). Because 
our noise has been propagated self-consistently, we use 
the noise model for our stacked sky-and-PSF-subtracted 
image to determine the la surface brightness limit for 
the ith aperture SBicr,i, and we take the mean value of 

l/N'^f SBicr,i of all apertures to be the surface bright- 
ness limit for the quasar pair SBi^r, which is quoted in 
Table 1. 
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Figure 3. Distribution of extracted surface brightness apertures 
relative to noise model predictions. The histogram shows the dis- 
tribution Xi = SBi /SBia^i for one of our sources SDSS J 0938+5317 
from Figure O In the absence of detectable extended Lya emission 
the SBi should be pure noise, and if our noise model is accurate, 
then Xi should follow a Gaussian distribution with unit variance, 
which is shown by the solid red curve. The standard deviation 
of Xi after clipping > 5cr outliers is a^hi = 0.97, which indicates 
that our noise model is indeed accurate and that fluctuations obey 
Gaussian statistics as expected. The median value of (Jchi, com- 
puted similarly, over all of our quasar pairs is 1.08, indicating that 
on average we tend to underestimate the noise by 8%. 

We can assess the accuracy of our noise model by 
examining the distribution of Xi = SBi/SBicr,i, where 
SB^ is the extracted flux of the ith aperture in each 
700 kms~^ X 1.0 arcsec^ window, and SBicr,i is the surface 
brightness limit based on our noise model at this loca- 
tion. In the absence of detectable extended Lya emission 
the SBi should be pure noise, and if our noise model is 
accurate, the Xi should thus obey a Gaussian distribu- 
tion with unit variance. The distribution of Xi foi" 
of our quasar pairs SDSSJ 0938+5317 is compared to a 
Gaussian in Figure [3l It is clear that our noise model 
provides an accurate description of the real noise in the 
stacked sky-and-PSF subtracted image. Systematic er- 
rors in the data reduction and sky and object modehng 
procedure are not accounted for in the purely statistical 
error SBicr,i, and thus we would expect to typically un- 
derestimate the true noise level. We can quantify these 
errors in our noise estimates by computing the standard 
deviation of Xi after clipping > 5cr outliers (resulting 
from e.g. a small number of residual cosmic ray hits or 
bad pixels that remain unmasked by our data reduction 
procedure). For the distribution in Figure [3] we obtain 
c^chi = 0.97, which indicates that for this case we actually 
slightly overestimated the noise. We similarly compute 
(Jchi for each quasar pair in our sample. The median value 
over all quasar pairs is 1.08, thus on average we tend to 
underestimate the noise by 8%. Because this is smaller 
than our estimated flux calibration error, we therefore 
quote statistical SB limits based only on our noise model 
as SBicr in Table 1. 

4.3. Visually Recovering Synthetic Sources 

Having quantified the formal errors in our sky-and- 
PSF-subtracted images, we now determine the signifi- 
cance level for a convincing detection. While we quote 
1(7 SB limits over 700kms~^ x l.Oarcsec^ apertures, 
our sensitivity to a given surface brightness is of course 
set by the size of the source and its spectral char- 
acteristics (i.e. kinematics). One can in principle al- 
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Figure 4. Illustration of detection significance of synthetic 
sources with various properties. The lower panel shows the Xskyonly 
image for SDSSJ 0938+5317 with synt hetic sources added accord- 
ing to the procedure described in ^4.31 In only the lower panel, the 
synthetic sources are noiseless, and their brightness and relative 
significance can be directly compared to the real (noisy) spectra of 
the two quasars shown in the central inset. The middle and upper 
panels show Xsky+object ctnd Xsmth images respectively, which are 
the same as the corresponding panels for SDSSJ 0938+5317 in Fig- 
ure [S] except that noisy synthetic sources have been added. The 
stretch and colormap are identical to that used for all of our x maps 
in Figure O and is indicated by the upper colorbar. The horizontal 
red and blue dashed lines and the vertical red dotted line are the 
same as in Figure O The colored circles indicate the location of 
the sources for various significance levels which are labeled in black 
at the top of the bottom panel. The red circles indicate the large 
'CGM' sources which are both spectrally and spatially resolved, 
black circles the 'LLS' sources which are spatially resolved but 
spectrally unresolved, and blue circles the 'CORE' sources which 
are both spatially and spectrally unresolved. This figure suggests 
that we should be able to marginally detect 'CGM' like sources 
that are > 2 x SBicr, 'LLS' like emission with > 3 x SBicr, and 
'CORE like sources to > 10 x SBio-. Note that the blob of emis- 
sion near Av = km s""*^ and impact parameter R± — 20 kpc is 
real small-scale Lyo; fuzz emission from the f/g quasar (see ^5.2|) . 



ways reach lower SB levels by averaging over larger 
spatial regions, and it is easier to detect narrow lines 
over the background. To address these issues we con- 
struct synthetic fluorescent sources and add them to 
the sky- and-PSF-subtr acted image for a representative 
pair SDSSJ 0938+5317. This pair has Zbg = 2.32, 
Zfg = 2.07, an angular separation of ^ = 5.6'^ cor- 
responding to R± = 49 kpc, and was observed with 
LRIS, as are the bulk of our pairs. The total exposure 
time was 1800s, allowing us to reach a surface bright- 
ness limit of SBicr = 0.29 x 10~^^ ergs~^ cm~^ arcsec"^ 
which is close to the median value of our sample SBi^r = 
0.31 X 10~^^ergs~^ cm~^ arcsec"^ (Table 1). Figure [3] 
testifies to the accuracy of our noise model for this object, 
and illustrates that the surface brightness noise fluctua- 
tions obey Gaussian statistics. Our spectral resolution 
for these observations is FWHM= 153 km s~^ and the 
seeing measured from the quasar profiles is 0.79^' corre- 
sponding to a spatial resolution of 7.0 kpc. 

We construct synthetic sources with a uniform spa- 
tial surface brightness distribution, where the size is 
specified by a single dimension along the slit (ifake- 
Spectrally, we assume a Gaussian velocity profile spec- 
ified by FWHMfake- The amplitude of the synthetic 
sources is set such that their total surface brightness 
(i.e. integrated out spectrally) is equal to SB = 
[1,2,3,4,5,10] X SBicr, which corresponds to SBLya = 



-200 [0.29,0.57,0.86,1.1,1.4,2.9] x 10"^^ ergs 



-1 —2 

^ cm ^ 



arcsec-2 for SDSSJ 0938+5317. 

We consider synthetic sources with three different spa- 
tial and kinematic distributions which represent our ex- 
pectations for the sources that might occur near a bright 
quasar. The first is a large (ifake = 75 kpc source with 
kinematics FWHMfake = 500 km s~^, which is both spa- 
tially and spectrally resolved. This source represents the 
kind of emission expected if the quasar CGM is made 
up of a large population of small spatially unresolved 
clouds which nevertheless have sufficient covering factor, 
/c, to result in detectable Lya emission (see Figure [TJd- 
d); we denote this model as 'CGM'. In this scenario, 
the kinematics of the emission correspond to the mo- 
tions traced by the metal-line absorbers as observed near 
quasars with b/g sightlines (Prochaska & Hennawi 2009, 
Prochaska et al., in prep.). The second synthetic source 
is meant to represent a spatially resolved Lya source 
(see Figure [1^), which is the counterpart to the fre- 
quently observed optically thick absorption in b/g quasar 
sightlines (Figure [2]). In this 'LLS' scenario, the size 
is taken to be (ifake = 20 kpc, corresponding to about 
three resolution elements. The kinematics are assumed 
to be comparable to the narrow turbulent and/or ther- 
mally broadened line widths ~ 20 km s~^ characteristic 
of metal absorption in individual LLS absorption com - 
plexes, (i.e., spectrally unresolved Prochter et al.|[2"010f ). 
For the third synthetic source, denoted as 'CORE', we 
assume a size (ifake of two LRIS pixels corresponding to 
4.8 kpc, and cold kinematics ~ 20km s~^, such that it 
is both spatially and spectrally unresolved. This source 
is meant to mimic the properties of the candidate flu- 
orescent emitters recently detected by Cantalupo et al.l 
(|2012l ) at ~Mpc distances from bright quasars. These 
fluorescent candidates tend to be compact, are spatially 
unresolved in ground-based seeing, and are believed to 
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be the dense cores of proto-galaxies. For spatially un- 
resolved sources, the total Lya luminosity is the rele- 
vant quantity rather than the surface brightness, and 
the corresponding luminosity of our CORE sources are 
LLya = [0.92,1.8,2.8,3.7,4.6,9.2] x 10"^^ ergs-^. Note 
that the median luminosity for the most compelling (i.e. 
those with W^ya > 240A) iCantalupo et al.l (j2Q12f ) fluo- 
rescent candidates is Li^ya = 9.1 x 10^^ ergs~^, which 
is at the level of the lOa model source we consider. Our 
model images are then convolved with the spectral and 
spatial resolution of our data, converted into counts us- 
ing our sensitivity function, and Poisson realizations of 
the average counts are added to the real sky-subtracted 
and sky-and-PSF-subtracted images. 

Figure H] illustrates the detection signiflcance of the 
synthetic source model profiles which have been added 
to the real data for SDSSJ 0938+5317. The lower 
panel shows the Xsky with synthetic sources added. 
In this panel only, the synthetic sources are noiseless, 
and their brightness and relative significance can be 
directly compared to the real (noisy) 2-d spectra of 
the two quasars shown in the central inset. Quantita- 
tively, a uniform screen of surface brightness at SBi^r = 
0.29 X 10~^^ erg s~^ cm~^ arcsec"^ emitted at the f/g 
quasar Lya with a top-hat velocity distribution matched 
to our 700 km s~^ spectral extraction aperture would re- 
sult in 0.44 e~ electrons per pixel, whereas the typical 
value of the la noise in our stacked images (far from the 
quasars) is a 4 e~ electrons per pixel. 

Inspection of Figure |4] suggests that we should be able 
to marginally detect 'CGM' like sources down to a level 
of > 2 X SBicr, 'LLS' like sources down to > 3 x SBicr, 
and the spatially and spectrally unresolved 'CORE' like 
sources to > 10 x SBicr. We emphasize that these sig- 
nificance limits are rough guidelines, as the noise level 
varies with wavelength due to fiuctuations in the sky 
level. For example, the reason why the 4(j sources in Fig- 
ure H] look scarcely more significant than the Scr ones, is 
that the 4cr sources happen to land on a region where the 
sky-background is ~ 10% higher. Furthermore, in some 
cases systematics in the data reduction limit our ability 
to make detections down to low significance levels. These 
systematics manifest as large-scale correlations in the x 
maps, and are particularly conspicuous in the smoothed 
Xsmth images. In the absence of emission, the Xsky+object 
image should be a random Poisson noise field because 
the pixel noise is completely uncorrelated, and thus the 
characteristic size of structures in the Xsmth should be 
comparable to the size of the smoothing kernel 5 — 6 
pixels. Some of our x maps in Figure [5] show correlated 
structure and spurious features due to systematics, the 
worst case being SDSSJ 0814+3250, where a large-scale 
gradient in the fiux level is seen across the spatial direc- 
tion. This large-scale gradient is likely due to a poorly 
determined slit-illumination function. 

5. RESULTS AND DISCUSSION 

In what follows we interpret our observational results 
from the previous section in terms of the emission mech- 
anisms discussed in ^ The results further constrain the 
physical properties of the CGM of the massive galax- 
ies hosting z ^ 2 quasars. Before delving into details, we 
briefiy summarize the main findings of our search for Lya 



emission in the regions surrounding quasars. Careful in- 
spection of the Lya emission maps in Figure [5] reveals 
the following: 

• Fluorescent Lya emission with properties expected 
in the 'mirror' approximation, is not detected 
at the locations of optically thick absorbers pin- 
pointed by the b/g quasar sight lines. 

• Despite the fact that the covering factor of optically 
thick absorbers on scales R± ~ 100 kpc is high 
(/c > 0.5), we do not detect diffuse Lya emission 
on these large scales, with the exception of a single 
system. 

• In a single system SDSSJ 0841+3921, we de- 
tect a spectacular Lya filament around the fore- 
ground quasar with surface brightness SB 

1 X 10~^^ergs~^ cm~^ arcsec"^ and an extent of 
230 kpc and possibly as large as 330 kpc(discussed 
further in Hennawi et al. in prep.). 

• Extended Lya fuzz is frequently detected on small- 
scales R± < 50 kpc. We detect definite fuzz in 10 
objects, or 34% of our sample, at levels of SBly^ — 
0.5 — 5 X 10~^^ergs~^ cm~^ arcsec"^. 

• In SDSSJ 0856+1158, a high-equivalent width 
{Wi^ya > 50A; rest-frame) Lya-emitter is detected 
with LLya =2.1 + 0.32 x 10^^ ergs" ^, at a distance 
of 134 kpc and velocity —540 km s~ relative to the 
f/g quasar. 

As we will frequently refer to our observations in ag- 
gregate, for reference we list the typical properties of our 
sample of 29 projected quasar pair observations. The 
median redshift of our f/g quasars is Zfg = 1.95, which 
are probed by b/g sightlines at a median impact pa- 
rameter of R± = 106 kpc. The luminosity at the Ly- 
man limit has a median value log^g ^^^ll — ^0.3 and 
our f/g quasars have a median de-reddened z-band ap- 
parent magnitude of 19.7. Our median la Lya surface 
brightness is SBi^r = 0.31 x 10~^^ergs~^ cm~^ arcsec"^, 
and according to the example in §4 3| large scale emis- 
sion from the quasar CGM should be detectable for > 
2 X SBi^ or SBLyc > 0.6 x 10"^^ erg s"^ cm-^ arcsec-^. 

5.1. Fluorescent Emission from Optically Thick Gas 

In this subsection we compare our observational con- 
straints on extended Lya emission to predictions for re- 
combination radiation from optically thick g as, as laid 
out in §2.2.21 We consider emission from spatially re- 
solved 'mirrors' and from an unresolved distribution of 
clouds in turn. 

5.1.1. Fluorescent Emission in the Mirror Approximation 

The frequent detection of optically th ick gas in 
the CGM of quasars quantified in §3.4[ an d dis- 
cussed i n other work in the QPQ seri es (Henna wi et al. 
^065 'H ennawi fc Prochaskal I2QQ7I : fProchaska et al. 
2012, 2013), provided initial optimism that we would de- 
tect significant fiuorescent Lya emission from this pop- 
ulation of clouds. In the simple 'mirror' approximation 
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(see ean.fn lG^ld fc WeinberilfToOGl) . ^ 60% of the ion- 
izing photons impinging on an opticahy thick cloud are 
converted into Lya photons, which may be 'reflected' 
into our Hne-of-sight. A precise estimate of the sur- 
face brightness depends on the ionizing photon flux, the 
cloud size, the geometrical configuration between the 
sourc e, observer, and cloud (s), as well as t he velocity 
field (jCantaluDO et al.ll2QQ5l : iKollmeier et al.ll2QQ8[ ). For 
spatially resolved clouds (see Figure [TJi), the Lya emis- 
sion should also be spatially resolved and emerge from 
the side of the absorber that lies closest to the quasar 
at a SB level given by eqn. ([12]). The predicted SBmirror 
values are given for each pair in Table 1, where the dis- 
tance is taken to be the impact parameter R± to the b/g 
si ght line. 

lAdelberger et all (j2QQ6f ) purported to detect Lya 
fluorescence from a serendipitously discovered DLA 
(logA^Hi = 20.4) situated 49" away from a lu- 
minous quasar at z = 2.84. The ionizing flux 
at transverse distance R± = 384 kpc implies a 
very high mirror surface brightness of SBmirror = 
1.1 X 10~^^ ergs~^ cm~^ arcsec"^ in this system, owing 
to the extremely luminous foreground quasar (r ~ 16). 
The authors detect a marginally resolved emission-line 
source near the b/g sightline with comparable surface 
brightness; and based on this agreement, the presence of 
a double-peaked kinematic emission line profile, the small 
spatial separation (1.5^^ or 11 kpc) of the emission from 
the DLA (and b/g sightline), and importantly, the fact 
that the emission is on the side of the absorber closest 
to the f/g quasar, they argue that it is Lya fluorescence 
from the same gas responsible for the DLA absorption. 

Here we have conducted a similar search for analogous 
'mirror' fluorescent emission from a much enlarged sam- 
ple of 29 si^htlin es with the same experimental setup as 
lAdelberger et al.l (j2QQ6f ). In 24/29 of t hese sigjhthnes we 
could have detected emission similar to lAdelberger etaTI 
(j2QQ6f ) at > lOcr (see Table 1), but fail to detect any- 
thing comparable. As discussed in §4.3| resolved sources 
with size o:^ 20 kpc, should be detectable down to the 
level corresponding to 2 x SBi^r. With the exception 
of SDSSJ 0841+3921, which we discuss further in Hen- 
nawi et al. (in prep.), we do not detect any emission in 
the immediate vicinity of the b/g quasar sight lines, as 
would be expected if the absorbing clouds were spatially 
resolved Lya mirrors, although we had the sensitivity 
to detect 24/29 sources at > lOa. Based on this large 
number of null detections of mirror flu orescence, we con- 
clude that the purported detection of lAdelberger et"al] 
()2006f ) is most likely the coincidence between a DLA 
and a star-forming galaxy, and was incorrectly inter- 
preted as Lya fluorescence. Indeed, the Lya emission 
in this source coincides with the location of a contin- 
uum source and the corresponding rest-frame equivalent 
width of only 75A is well within the range produced by 
conventional stellar populations (i.e. < 240 A), and is 
hence entirely consistent with being powered by star- 
formation. It is possible that the Lya emission in this 
object is boosted by a fluorescent contrib ution, an effect 
also discussed in iCantalupo et al.l ()2012f ) , but given its 
relatively low rest-frame equivalent width, and our large 
sa mple of similar nul l dete ctions, any interpretation of 
the lAdelberger et al.l ()2006l) object as Lya fluorescence 



is entirely speculative. 

5.1.2. Fluorescent Emission from Unresolved Clouds 

For unresolved optically thick clouds distributed 
around the quasar (see Figure [Tb-d), the fluorescent 
emission will be manifest as an extended and spatially 
resolved 'fuzz', similar to the large Lya nebulae observed 
in HzRGs and Lya blobs. In this case, the average sur- 
face brightness is given by eqn. ([15]). Given the high 
fc ^ 0.50 covering factor of optically thick gas at dis - 
tances < 200 kpc (see Figure [2 Prochaska et al.] l2012f ). 
the expected surface brightness of such nebulae is about 

30% smaller than the mirror values tabulated in Ta- 
ble 1 averaged over an aperture R±. However a RJ^ 
profile is expected for the idealized case of a uniform 
distribution in a spherical halo (see eqn. [15]), and hence 
the emission would be brighter closer to the quasar. We 
emphasize that if the optically thick gas giving rise to 
the high covering factor is illuminated by the f/g quasar, 
the expectation is that a SB comparable to the mir- 
ror sufrace brightness is expected at impact parameters 
R± < 200 kpc all along the slit (see Figure [J). 

Such large-scale diffuse Lya halos would have been 
very easy to detect. Indeed our CGM synthetic source 
simulation with emission extending over 70 kpc (see § 14.31 
and Figure |4]) indicates that we could have detected such 
resolved emission at levels SBLya ^ 2 x SBi^r. Assum- 
ing a fixed distance R = 100 kpc and a covering fac- 
tor of fc = 0.50, consistent with our measurements in 
Figure [21 the optically thick surface brightness given by 
eqn. ([T5]) is > 5 x SBi^r for 28/29 targets in our sample, 
and > 10 X SBi^r for 19/29 targets. No diffuse Lya emis- 
sion is detected along the slit at distances of > 50 kpc 
from any of our f/g quasar, with the exception of the 
large-scale filament in SDSSJ 0841+3921, which we dis- 
cuss further in Hennawi et al. (in prep.). 

5.1.3. Why is Optically Thick Fluorescence Absent? 

In the preceding subsections we argued that our data 
shows no evidence for fluorescent emission from optically 
thick clouds. The possible explanations for these null 
results include: 

• The optically thick clouds have very small radii. 

• The gas lies at much greater distance from the 
quasar than implied by the observed impact pa- 
rameters. 

• The emission is extincted by dust grains. 

• The optically thick gas is not illuminated by the 
quasar. 

The first point about the clouds sizes relates to the 
fact that in the mirror approximation (eqn. [T2]) the Lya 
flux at the location of the b/g sightline is proportional to 
the absorber cross-sectional area ctc, fhya ^ SBLya x (Jc. 
As noted above, we would have easily detected emission 
from resolved sources with cloud sizes > 10 kpc. But if 
the clouds are very small, such individual mirrors will not 
be detectable. Indeed, our detailed absorption line mod- 
eling of one member of our sample SDSSJ 1204+0221 in 
QPQ3, yielded very small cloud sizes ~ 10 — 100 pc, rel- 
ative to the typical ~ 5 kpc resolution implied by our 
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ground based seeing, although it is stih unclear whether 
this one system represents the typical properties of gas 
in the quasar CGM. Nevertheless, even for such unre- 
solved clouds (cloud sizes <C 10 kpc) the high covering 
factor, /c, of optically thick gas still implies significant 
resolved Lya fiux, resulting from the collective emission 
of a population of unresolved clouds (see Figure[T]). Thus, 
while an appeal to small clouds can explain our null de- 
tections of Lya mirror fiuorescence from the individual 
cloud absorbing the background sightline, a large SB is 
still expected from unresolved clouds. 

The second explanation is that we have substantially 
over-estimated SBmirror because we have assumed the gas 
is at distance R± (eqn. [T2]) . whereas in reality it lies 

at a significantly larger distance r = + j_ , where 

R\\ is the unknown line-of-sight distance. This would 
reduce our estimates for the mirror SB by [1 + {R\\/ 
R±)'^]^ a large factor if R\\ is significantly larger than 
R±. We consider this interpretation highly unlikely 
for several reasons. First, note that we have demon- 
strated that the covering factor of optically thick ab- 
sorbers is hi^h at imp act paramet ers R± < 200 kpc from 
quasars (iHennawi et al. 2006ijH ennawi fc Prochaskal 
120071 : iProchaska eTal.l .2012. I2013L see e.g. ^Ti]) . but 
we observe a decrease for R± > 200 kpc, both for our 
parent sample shown in Figure [2] as well as in larger 
samples (Prochaska et al. 2012, 2013). The presence of 
this dropoff at R± 200 kpc indicates that we are in- 
deed probing gas interior to the quasar halo on scales 
r < 200 kpc. Second, our absorption line analysis of one 
member of our sample in IProchaska fc Hennaw i' (2009) 
revealed extreme kinematics and a very high metallicity 
in comparison to the typical properties of LLSs, imply- 
ing that this gas is local to the very atypical f/g quasar 
environment. We also find large equivalent widths for 
low- ion transitions (e.g. C II 1334) for the parent sam- 
ple (jProchaska et al.ir2012f ). which are uncharacteristic of 
optically thick absorbers along random quasar sight lines. 

It is conceivable however, that optically thick absorbers 
very close to the quasar are photoionized away, resulting 
in a donut-like radial distribution. In this case, the ab- 
sorbers in our sample with the the smallest impact pa- 
rameters R± < 50 kpc and highest corresponding mir- 
ror surface brightness, could actually be at much larger 
physical distances of ~ 200 kpc, reducing the anticipated 
SB by factors > 10. A clustering analysis of our full 
sample of optical ly thick absorbers aro und quasars will 
be presented in (Prochaska et al.|[2013L see also QPQ2); 
there we explore the degree to which donut like distribu- 
tions can be ruled out. Note, however, that we do ob- 
serve evidence for an increased frequency of the highest 
column density absorbers (sub-DLAs and DLAs) at the 
smallest impact param eters R < 50 kpc, (QPQ1,QPQ2 
IProchaska et al.|[2013f ): the obvious interpretation is that 
this results from increased covering factor and/or col- 
umn density at smaller radii, again implying that the 
true radius is comparable to the impact parameter. Fi- 
nally, note that even if a donut like radial distribution 
implies that we have overestimated the SBs for the small- 
est impact parameters {R < 50 kpc), we still could 
have detected Lya emission from optically thick gas at 
R = 100 kpc in 28 out of 29 members of our sample 



at high statistical significance. Based on all of these ar- 
guments, we conclude that radial distances much larger 
than the impact parameters are unlikely, and therefore 
do not explain our null detections. 

In the Appendix we consider the third explanation, 
that resonantly trapped Lya photons are extincted by 
dust grains, in detail. There we show that, given 
reasonable assumptions about the column density and 
ionization state of absorbers in the quasar CGM, and 
rather conservative assumptions about the metallicity 
and hence dust content of this gas, extinction by dust 
grains can result in a factor of at most ~ 2 reduction 
in the total Lya emission. Furthermore, this estimate 
applies to the case of Lya photons generated uniformly 
throughout the volume of an emitting cloud, as it would 
be for collisionally-excited cooling radiation. It thus 
likely overestimates the attenuation of fluorescent Lya 
photons from an externally illuminated optically thick 
absorber. In these self-shielding clouds, Lya photons 
traverse a smaller optical depth before escaping than we 
assumed in the Appendix. 

Having cast significant doubt on the other three in- 
terpretations, we consider the last possibility, that the 
optically thick gas detected in b/g sightlines is not illu- 
minated by the quasar, as the most probable explanation 
for our null detections. In principle this could result from 
either of two physical effects, anisotropic emission of the 
quasars ionizing radiation, or intermittent emission by 
the quasar on timescales comparable to the light cross- 
ing time tcross = R±/c ~ 3 X 10^ yr for R± ~ 100 kpc 
(see also discussion in QPQ2). Note however that the 
frequent detection of extended Lya fuzz at distances 
R± ~ 30 kpc (see §5.2p indicates that quasars emit con- 
tinuously for at least Across = 30 kpc/c, hence under the 
intermittency interpretation the duration of quasar emis- 
sion bursts would have to follow a very narrow distribu- 
tion, which seems rather unphysical. 

We thus favor the anisotropic emission interpretation, 
which results if, e.g. the f/g quasar is obscured from 
some vantage points, and the optically thick gas we do 
detect in b/g sightlines is shadowed. In addition to the 
absence of fluorescent Lya emission, several independent 
pieces of evidence support this hypothes is: (1) in the con- 
text of uni fled models of AGN (see e.g. lAntonuc"cilll993l : 
lElvisll2QQQh . there is a wide and growing body of evidence 
for a signiflcant population of obscured quasars bot h in 
the local and distant Universe (e.g. Willott et al . 20001: 
Treister & Urrvl 120051: iGilh et a l. 2007; Maioli no et al l 

.2007: .Treister et al.ll2008l : iReves et al 2008: Stern et al.l 

l2012f ). which implies quasars emit their UV radiation 
anisotropically; (2) our comparison of the incidence of 
optically thick gas transverse to and along the line-of- 
sight to quasars indicates an anisotropic clustering pat- 
tern (see Figures gMl and also QPQ2); (3) the phys- 
ical conditions of the absorbing gas in one of our sys- 
tems SDSSJ 1204+0221 determined from detailed ab- 
sorption line modeling, does not show strong evidence 
that the f/g quasar is shining on it (QPQ3); (4) we flnd 
a large ratio of low-ion to high-ion equivalent widths 
(C II/ C IV) for the f/g quasars in Table 1 as well as 
in the parent sample (QPQ5). These values are com- 
parable to typical ratios in DLAs, suggesting a low- 
ionization state; (5) at distances of r ^ 100 kpc, sev- 
eral studies indicate that the ionizing radiation emis- 
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sion from the qu asar will 
est gas , 



tne qu asar will photo i onize all but the dens- 
Cantal upo et alJ I2QQ5I : iHennawi fc Prochaskal 
2QQ7C iCh elouche et aL' '2QQ7 I: iKollmeier et all I2QQ8I : 
Hennawi et al. 2009; CantaluD oet al.ll2012[ ). 



If the quasar emits its radiation anisotropically, we can 
crudely estimate the solid angle of the emission ft, and 
consequently the implied fraction of AGN which are ob- 
scured /obscured = 1 — ^/47r. In the picture we have 
put forth, clouds detected as optically thick absorbers 
are not illuminated and hence do not emit Lya, whereas 
clouds which are illuminated are photoionized, optically 
thin, and although they emit Lya radiation, it is n ot de - 
tectable at our current sensitivity (see discussion in §5.2p . 
Note that in this context, the covering factor that we es- 
timated in §3.41 is actually the intrinsic covering factor 
of optically thick clouds, /c, reduced by the impact of 
photoionization. If we assume that the quasar emits into 
a solid angle then in our simple model of a uniform 
distribution of clouds (see §2.ip . we can rewrite eqn. ([1]) 
for the observed covering factor /c,obs = /c(l — r^/47r), 
where fc now represents the intrinsic covering factor, 
e.g. before the quasar turns on. This expression assumes 
that all illuminated clouds are photoionized and is valid 
provided that the intrinsic fc < lE3- As we do not have 
knowledge of the intrinsic covering factor fc , we can only 
obtain an upper limit on the opening angle, which as- 
sumes that all of the 'misses' in Figure [2] occur because 
our b/g sightline intersects the volume which has been 
hollowed out by photoionization. This limit corresponds 
to setting the intrinsic covering factor fc to unity. Then 
if we take /c,obs = 0.5 consistent with our measurements 
(Figure [2j), we obtain 0/47r < 0.5. The obscured fraction 
of z ^ 2 quasars is then constrained to be /obscured > 0.5 
which is broadly consistent with independent estimates 
of the obscured fraction from mult i- wavelength st udies 
of AGN (e.g. iWillott et al . 2000; Trei ster &: Urrvl [2005 



Gilli et al.l[2007[ ? :lTreister et al 200& iReves et al.ll2008 ; 

Stern et al.l 120121 : iLusso et al.l 120131 ). However, we em- 



phasize that our model for cloud photoionization is ex- 
tremely crude, and surveys for AGN at different wave- 
lengths (X-ray, optical, near/mid-IR, radio) often come 
to different conclu sions about obscured fractions (see e.g. 
ILusso et al.1l20l3[ ). 

5.2. Small-Scale R± < 50 kpc Lya Fuzz 

In the previous section, we argued that the absence 
of significant Lya emission from optically thick gas in 
the f/g quasar CGM indicates that this material is not 
illuminated by the quasar. However, our longslit obser- 
vations frequently reveal extended Lya fuzz at a level of 
SB ~ 0.5 — 5 X 10~^^ergs~^ cm~^ arcsec"^ at distances 
R± ^ 50 kpc. To quantify the frequency of this extended 
emission, both authors independently visually inspected 
the maps in Figure [5] to search for compelling extended 
Lya emission. Each object in our sample (see Table 1) 
was classified as having definite fuzz, ambiguous, or hav- 
ing no fuzz. Classifications were merged according to the 
following rules: 1) an object is classified as definite or no 
fuzz only if both authors agreed on that classification; 
2) all other objects were considered ambiguous. Out of 

12 por values of fc significantly larger than unity, this expression 
breaks down because one needs to take into account the constraint 
that the quasar is always pointing towards the observer. 



our sample of 29 systems, this yields 10 objects with def- 
inite fuzz, 8 ambiguous objects, and 11 for which no fuzz 
is detected. Given that this fuzz tends to be asymmet- 
rically distributed, and that our single-slit orientation 
covers only a fraction of the area around the quasar, the 
true detection rate of such fuzz is expected to be sig- 
nificantly higher. Our detection rate is hence about a 
factor of two lower than the frequency of ~ 50 — 70% 
deduced fr om the more recen t studies of this small-scale 
Lya fu zz (Christensen et al. 



2006; Courbin et al.l 120081 : 
IHenna wi et al. 2009; North et al. 2012), although lim- 
ited statistics, sample inhomogenities, differing method- 
ologies, and our single position angle preclude a quanti- 
tative comparison. 

We now interpret this small-scale Lya fuzz in terms of 
the emission mechanisms discussed in ^ Consider our 
typical f/g quasar at z 2, which has log^^g ^z^ll 

0^ 30.3. 

If we assume R_\_ = 35 kpc as characteristic of the 
impact parameters at which we detect the fuzz, the 
average surface brightness for optically thick fiuores- 
cence given by eqn. ([15]) is SBLya = 1.0 x 10~^^(/c/ 
1.0) ergs~^ cm~^ arcsec"^. It is thus quite clear that 
small-scale fuzz cannot be arising from the same pop- 
ulation of clouds giving rise to the high covering factor 
/c ^ 0.50 of optically thick absorbers at R± < 50 kpc in 
Figure [21 since the implied surface brightness would be 
two orders of magnitude higher. To be consistent with 
optically thick emission, the covering factor of clouds re- 
sponsible for the small scale fuzz needs to instead be 
much lower fc — 0.005 — 0.05, and this low covering 
factor matches the values deduced in prev ious work on 
quasar nebulae (e.g. iHeckman et al.lll991bl ). 

Considering together both the detected small scale fuzz 
emission, which implies a small covering factor fc — 
0.005 — 0.05, and the lack of large-scale nebulosity from 
gas with a high fc ^ 0.50 covering factor established by 
b/g sight lines, the picture that emerges is that clouds 
with a range of densities and covering factors exist in 
the CGM around quasars. In the absence of ionizing 
photons, a low density nn ~ 0.01 — 1 cm~^ population 
dominates the covering factor, and it is this gas that 
we typically observe in absorption with b/g sight lines. 
In contrast, the denser clouds (nn ~ 10 — 100 cm~^) 
cover only a tiny fraction of the line-of-sight ~ 0.01. 
When illuminated by the f/g quasar, the lower den- 
sity (high covering factor) clouds are photoionized by 
the intense quasar radiation, whereas the dense clouds 
(low covering factor) can self-shield and survive. It 
is then these dense clouds which are responsible for 
the small scale Lya fuzz that we detect in our sam- 
ple, and whic h have been prev iously detected a round 



quasars (e.^. Bunker et aTl 120031: IWeidinger et al] 



Christensen et al.l I2Q06I; [F rancis & McDonnell 

Courbin et all 120081 : IHennawi et al. 2009; Nor th et all 
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2012[ ). This scenario could also apply to the Lya 



blobs and HzRGs. Indeed, it is intriguing that detailed 
modeling of the spectra of extended emission-line re- 
gions (EELRs) around quasars and FR II radio galax- 
ies at low redshift z < 0.5, similarly show evidence for 
two phase media composed of an abundant low den- 
sity medium (nn ~ 1 cm~^) and much rarer high den- 
sity clouds (ng ^ 400 cmj^ iStockton et al .'' 2002, 20^ 
iFu Stockt^ I2006L I2007L l2009f ). It thus seems plan- 
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sible that the smah-scale fuzz that we detect in Lya 
are high-redshift analogs of the E ELRs detected around 
0.5 type- 1 (e.g. iStockton et al. 



low redshift z < ^.^ , _ 
iHusemann et al.|[2Q12f ) and tvpe-II (fCreene et al 



200' 



20111 



quasars, in [O III] and Balmer lines. 

While the explanation of small-scale fuzz above is im- 
minently plausible, and has indeed been the standard 
interpretation for nebulosities around both HzRGs (see 
e.g. [McCarth y et aLlpOOl: lMcCarthvllf993l) and quasars 
(e.g. Heckma n et al. 1991b), as emphasized in ^2.4[ the 
interpretation of Lya emission alone is degenerate, and 
optically thin fluorescence and Lya scattering can also 
result in comparable emission levels. First consider opti- 
cally thin fluorescence from a population of clouds with 
(i? , nn , A^H , /c) = (35 kpc , 1.0 cm-^ , lO^^-^ cm-^ , 1.0). 
This results in a SBLya = 10~^^ergs~^ cm~^ arcsec"^ 
according to eqn. (p!Q|) . and it is rather intriguing that this 
column density Nn and volume density nn are compara- 
ble to values deduced from our detailed absorption line 
modeling of SDSSJ 1204+0221 QPQ3, and the covering 
factor consistent with the large values we deduced in §3.41 
In other words, if such clouds are not illuminated by the 
quasar, they could account for the optically thick absorp- 
tion that we observe in b/g sightlines at R± < 200 kpc. 
But if illuminated at a distance R = 3b kpc, such clouds 
would be so highly ionized that they would appear op- 
tically thin in absorption (e.g. along the line-of-sight). 
Hence, the small scale fuzz could be arising from the il- 
luminated (hence optically thin) analogs of the shadowed 
(hence optically thick) clouds that we detect in absorp- 
tion in background sightlines. 

But if this interpretation is correct, why don't we see 
fuzz at larger distances R± 100 kpc? Recall that in 
the optically thin regime, fluorescence does not depend 
on distance provided that the clouds remain optically 
thin, and at R± :^ 100 kpc this would marginally be the 
case (logA^Hi — 17.2). Thus under this interpretation, 
there should be extended emission at R_\_ 100 kpc 
which is just as bright as the small-scale fuzz. Such 
emission is definitely not observed in our dataset. As 
the optically thin surface brightness scales with the com- 
bination /cA^H^H, one solution would be to make this 
product substantially lower at ~ 100 kpc compared 
to R ^ 35 kpc, so that the optically thin emission would 
then be too faint to be detected. 

Now consider the SB levels from Lya scat- 
tering. A model with (i^ , , A/'h , /c) = 
(35 kpc,0.1cm-^,10^°-^cm-2,1.0), i.e. nearly identi- 
cal to the optically thin case considered above, but with 
an order of magnitude lower volume density nn, results 
in scattering SBLya = 0.9 x 10~^^ergs~^ cm~^ arcsec"^ 
(see eqn. [20]) which is also bright enough to be the mech- 
anism responsible for the small-scale Lyo; fuzz. Further- 
more, the lower volume density nn of this configuration 
lowers the optically thin fluorescence discussed above 
to SBLya = 8 X 10~^^ ergs~^ cm~^ arcsec"^, almost an 
order of magnitude fainter than our typical 2a detection 
limit SBLya 6 X 10~^^ ergs~^ cm~^ arcsec"^. 

Thus both optically thick and optically thin fluores- 
cence, as well as Lya scattering can all plausibly explain 
the SB levels of the small-scale fuzz observed at distances 
R± < 50 kpc from quasars. Our simple estimates are 
crude, and more detailed numerical treatments which 



properly treat ionizing and resonant radiative transfer 
are required to improve their accuracy as well as make 
more detailed predictions about the SB profiles and the 
kinematics of the emission from these various processes. 
However, we believe that the key point, that the interpre- 
tation of the Lya emission alone is not unique, is robust. 
We have argued that additional information about the 
gas distribution from absorption lines measurements can 
ameliorate this issue. For example, if the covering fac- 
tor, column densities, and volume densities of the gas 
at < 50 kpc can be statistically characterized, as we 
have undertaken in the QPQ series, then we can make 
better predictions for the optically thin and scattering 
surface brightness discussed above. Moreover, observa- 
tions of additional emission line diagnostics, such as a 
non-resonant Balmer line (i.e. Ha), could distinguish 
between fluorescence and Lya resonant scattering. The 
possibility will always remain that a very low covering 
factor component, which cannot be easily mapped with 
absorption lines, dominates the emission, as is the case 
for the low covering factor optically thick fluorescence 
discussed at the beginning of this section. Nevertheless, 
a comparison of absorption line properties to emission 
line measurements yields important constraints on the 
physical properties of the gas dominating the covering 
factor. 

To conclude our discussion of the small-scale fuzz, we 
note that one can in principle also set a lower limit on the 
opening angle of the quasar emission based on the fact 
that fuzz is frequently detected extending to a distance 
R^ 50 kpc. Unfortunately, given the unknown dis- 
tance r to the gas, this is not particularly constraining. If 
we conservatively assume that the fuzz-emitting gas has 
to lie within the quasar halo at r < 200 kpc, then sim- 
ple geometric considerations imply an opening angle Q/ 
47r > 0.004 (/obscured < 0.996) corresponding to emission 
into a cone with half-angle > 7°. If we instead assume 
the unknown line-of-sight distance R\\ is comparable to 
R±^ then we obtain Q/Att > 0.04 (/obscured < 0.96) cor- 
responding to a cone with half-angle > 23°. 

5.3. Comparison to Other Searches for Lya 
Fluorescence: The Lya Emitter in 
SDSSJ 0856^1158 

5.3.1. Previous Work on Lya Fluorescence 

[Cantalupo et al.[ (j2012L see also Cantalupo et al. 2007) 
conducted a survey for fluorescent Lya emission via a 
narrow band imaging survey of a 3.5 Mpc x 3.5 Mpc field 



around a very luminou s quasar with log^^g 



31.64 



(|D'Odorico et al.ll2008[ ). which is 20 times brighter than 
the typical f/g quasar in our sample. They detect an 
enhancement of sources with high rest-frame equivalent 
width EWLya > 240 A around the quasar, as compared 
to the 'field', which they interpret as the signature of Lya 
fluorescence. The significance of the 240 A rest-frame 
EW limit is that normal stellar populations are not ex- 
pected to produce EWs this large. These fluorescent can- 
didates tend to be compa ct and are spatia ll y unr esolved 
in ground based seeing. [Cantalupo et al.[ (|2012[ ) argue 
that they are the dense cores of proto-galaxies, which are 
either sufficiently dense to self-shield the quasar ionizing 
radiation and act as mirrors, or which are highly- ionized 
but sufficiently dense to emit observable optically thin 
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Lya. 

Given the different search strategy (narrow band imag- 
ing versus sht spectroscopy), the fact that their obser- 
vations reach a rest- frame SB0 a factor of 2.5 deeper 
than our median hmit, the much brighter quasar that 
was observed, and the much larger field-of-view they sur- 
veyed, a detailed com parison of our observations to the 
iCantalupo e t al.l ()2Q12i) results is not straightforward. 
The biggest challenge is that the vast majority of the fluo - 
rescent candidates discovered by 'Cant alupo et al.l (|2012') 
are at Mpc distances corresponding to several arcminute 
separations; whereas, our longslit maps in Figure [5] typ- 
ically extend only to an impact para meter of R± :^ 
200 k pc or about 25'^ However, note ICantalupo et al.l 
(|2012f ) did detect one source very near the quasar with 
= 15" corresponding to R_\_ = 127 kpc, which overlaps 
with the scales covered by our observations. This source 
had the highest rest-frame Wi^ya > 483 A in their sur- 
vey, and its luminosity Li^ya = 9.1 ± 0.7 x lO^^ergs"^ 
puts it among a handful of the brightest sources discov- 
ered, making it one of the most compelling fluorescent 
candidates. Consider now the total effective area cov- 
ered by our aggregate slit-spectroscopic observations of 
29 f/g quasars: Aeff = 29 x 150 kpc x 8.7 kpc. Here we 
compute the total area interior to a radius of 150 kpc 
(comparable to the distance of the Cantalupo et al. flu- 
orescent emitter) and covered by the 8.7 kpc width of 
our 1.0^' slits {z = 2). This area is 1.07 times the to- 
tal area 7r(150 kpc)^ around a single quasar. Hence we 
are motivated to search our emission maps around fore- 
ground quasars for compact fluorescing sources near the 
f/g quasar, similar to that discovered bv ICantalupo et al.l 
([2012). 

5.3.2. The Lya Emitter in SDSSJ 0856+1158 

Careful visual inspection of the emission maps in 
Figure [5] reveals statistically significant Lya emission 
from a compact object at a distance of —15^' from the 
f/g quasar in SDSSJ 0856+1158, which translates to 
= 131 kpc at the redshift of the Lya emission line 
z = 1.759. Extracting a spectral aperture at the Lya 
emitter (LAE) location with FWHM set to the seeing 
of the observations, we measure the LAE Lya flux to 
be /Lya = 0.94 ±0.14 x 10~^^erg s~^ cm~^, correspond- 
ing to Li^ya = 2.1 ± 0.32 X lO^^ergs"^. A gaussian 
fit to the extracted flux gives a velocity dispersion of 
ay = 147 km s~^, which corresponds to :^ 2.2 spectral 
resolution elements and is hence marginally spectrally 
resolved. A similar analysis of the spatial profile gives 
a FWHM within 20 - 50% of the 0.97'' seeing, which 
corresponds to 8.5 kpc at the redshift of Lya. Given 
the low signal-to- noise ratio, we assume the source is 
spatially unresolved, although it could be marginally re- 
solved. The velocity difference between the LAE and the 
f/g quasar redshift is 870km s~"^, but this large shift may 
also have a contribution from uncertainties in the quasar 
redshift, which is derived from the Mg II emission line 
and has an associated uncertainty of <Jz = 272 km s~^ 

When comparing the ICantalupo et al.l (|2012l ) results at z = 
2.4 to our results at z 2.0 we refer to the rest-frame SB 
(ergs"-*^ cm~^kpc~^) which scales as SBrest oc SBobs(l + z)'^/Da 
where SBobs is the observed SB and Da is the angular diameter 
distance. 



([Richards et al.|[2Q02b[ ). The velocity difference between 
the LAE and the metal lines in the DLA (coincident with 
the f/g quasar) is —537 km s~^, which is consistent with 
the kinematics expected in the massive dark matter ha- 
los hosting quasars (see e.g. QPQ3) as well as the typical 
metal-line kinematics seen in the quasar CGM using b/g 
sighthnes (Prochaska et al., in prep.). 

No continuum emission is detectable by eye at the 
location of the LAE over the entire spectrum. To ob- 
tain quantitative constraints on the continuum level and 
hence the rest-frame EW, we consider a 200 A region 
at the red end of LRIS-B spectrum (3668-3868A) where 
the observations are most sensitive. The average sin- 
gle pixel la error of our extracted spectrum is ax = 

4x 10~^^ ergs~^ cm~^ A , corresponding to an AB 
magnitude mareg = 25.7, where 3769 A is the mean wave- 
length over the region considered. Formally, the error on 
the mean continuum averaged down over the npix = 471 
pixels in this 200 A stretch is acont — cta/v^^Vx- This 

value is (Jcont = 2.0 x 10~^^ ergs~^ cm~^ A ^ or an AB 
magnitude of msred = 29.0, and the corresponding la 
lower limit on the rest-frame equivalent width would be 
Wi^ya ^ /Lya/c^cont > 173 A. In practicc, systematic er- 
rors in the spectroscopic data reduction probably prevent 
us from achieving this formal sensitivity to the mean con- 
tinuum, and hence we would not actually have been able 
to detect a 29th magnitude source ^ 20 times fainter 
than the ax flux level. Note however that the mean ex- 
tracted flux over this spectral region at the location of the 
LAE is 0.04crcont, and this fluctuation is at least consis- 
tent with being drawn from a Gaussian distribution cen- 
tered on zero with standard deviation (Jcont, i-e. the data 
does not show strong evidence for systematic error in the 
average flux level. We nevertheless adopt W^ya > 50A as 
a conservative lower limit on rest-frame equivalent width, 
which corresponds to using roughly 3.5crcont as the con- 
tinuum value, or assuming that we could have detected 
a source four times fainter the than ax pixel noise corre- 
sponding to m3769 = 27.6. 

How likely is the discovery of a high Wi^ya LAE within 
15^' and at small velocity separation from a quasar at 2: ~ 
2? The answer clearly depends on the physical process 
powering the Lya and we consider two possibilities: 

• The Lya is powered by star- formation just as for 
the general LAE population. Clustering enhances 
the probability of finding such galaxies close to a 
quasar. 

• The Lya emission is fluorescence powered by the 
f/g quasar's ionizing radiation. 

We consider both possibilities in turn. 

5.3.3. Clustering of LAEs around Quasars 

The expected number of LAEs with rest-frame EW > 
Wi^ya within a volume V around a quasar can be written 

A^LAe(> W^Lya) = nLAEP(> W^Lya) / dV [1 + ^QG{r)] , 

Jv 

(23) 

where tilae is the average number density above the 
canonical rest-frame EW threshold Wi^ya > 20A for se- 
lecting LAEs, P(> Wi^ya) is the cumulative distribution 
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of LAE EWs, and the integration over the correlation 
function ^qq quantifies the increase due to quasar-LAE 
clustering. 

Although there is strong evidence that the power- 
ful HzRGs trace large proto-clust er like ov erdensities at 
z>2 (e.g.Pen t ericc i et al. 2000: E) e Breuck et al. 2004; 
IVenemans eral.ll2005l : iKuiper et al.ll2012[ ). the clustering; 
of galaxies around more typical radio-loud quasars like 
SDSSJ 0856+1158 and/or the radio-quiet quasars which 
dominate our sample, appears to be much more mod- 
est (e.g. Fvnbo et al. 2001; Adelberger & Steidel 2005; 
fXrainor fc: Steidell 1201 2). We thus assume that the clus- 
tering of LAEs around quasars at 2; ~ 2 is the same 
as that of t he continuum sele ct ed LB G population at 
z - 2.7. ' Trainor fc Steidell (j2012[ ) (see also Adel- 
berger & Steidel 2005) measured the cross-correlation 
between LBGs and quasars at z ~ 2.7, and assuming 
a power law form ^ = (r/ro,QG)~^ with 7 = 1.5, they 
measured a quasar-galaxy cross-correlatiorQ length of 
^o,QG = 7.3/1""*^ Mpc (comoving). We integrate this 
correlation function over a cylindrical volume defined 
by an annular region [i^±,min, ^±,max] and a line-of- 
sight extent /[H{z)\. The minimum radius is set 
to i^±,min = 50 kpc (proper), motivated by the fact 
that it could be difficult to distinguish LAE emitters 
from the ubiquitous small-scale Lya fuzz seen around 
quasars. The maximum radius is set to i^±,max = 
130 kpc (proper; 15'^ at z 2), which is approximately 
the impact parameter of the LAE in SDSSJ 0856+1158 
and the bri ght ffuorescence c a ndida te discovered near the 
quasar by iCantalupo et al.l ()2012l ). Along the line-of- 
sight, we take Av = 1500 km s~ , which encompasses 
both gravitational motions of galaxies in the quasar halo, 
quasar redshift errors, and also agre es with the inter- 
val us ed in the clustering analysis of iTrainor fc SteideS 
(|2012|). Furthermor e , the 4nm narrow-band filter used 
by 'Cantalup o et al.l ()2012f ) also translates to a line-of- 
sight velocity interval of 2900 km s~ or about 2 x Av. 
Based on these assumptions, we find that over this vol- 
ume, clustering alone enhances the number density of 
LAEs around quasars by a factor of 19 above the cosmic 
average. 

For the background number density of LAEs nLAE, 
we first need to decide which Lya luminosity thresh- 
old to consider. The LAE in SDSSJ 0856+1158 has 
^Lya = 2.1 X lO^^ergs"^, but owing to the long- 
exposure (4800s) and relatively low redshift z = 1.76, 
this represents the most sensitive observation in our 
sample. As we wish to consider our sample in ag- 
gregate, we scale this value up by a factor of three, 
which represents the ratio of the luminosity limit in 
SDSSJ 0856+1158 to the median value of our sample. 
This value of Li^ya = 6.4 x lO^^ergs"^ is compara- 
ble to the luminosity of the fiuorescent candidate de- 
tected closest to the q uasar Li^y^ = 9.1 x lO^^ergs"^ 
from iCantalupo et al.l (j2012l ) , and also matches the me- 
dian luminosity of their best fiuorescence candidates. 

1^ Although the Trainor & Steidel (2012) study focused on hy- 
perluminous quasars, whereas the f/g quasar we consider here are 
not hyper luminous, there is no compelling evidence that z ~ 2 — 3 
quasar clustering strength varies with luminosity in either the auto- 
correlation (Shen et al. 2008) or the cross-correlation with galaxies 
(tTrainor &; Steidel 2012; Adelberger & Steidel 2005). 



We use the iCiardullo et al.l ()2012l ) Schechter luminos- 
ity function fi ts to the z = 2.1 LAE survey data of 
iGuaita et al.l ()2011f ). to determine that the comoving 
number density of LAEs with Li^y^ > 6.4 x 10^^ ergs~^ 
is nLAE = 5.1 X 10-^/i^Mpc-^ ICiardullo et^ (j20T2h 
also estimated the rest-frame equivalent width distribu- 
tion, and found that z = 2.1 LAEs selected to have 
rest-frame T^Lya > 20A, follow an exponential distribu- 
tion with rest-frame scale length Wq = 50A. Hence the 
cumulative probability distribution of rest-frame EW is 
^(> W^ya) = exp [-{W^ya - 20A)/VI^o]. Given our con- 
servative lower limit for the rest-frame WLya > 50A of 
the LAE in SDSSJ 0856+1158 we obtain a probability 
P(> 50 A) = 0.55, whereas our formal limit of Wj^ya > 
173A gives P(> 173 A) = 0.047. Putting ah the pieces 
together (see eqn. l23|) . we find that the expected number 
of LAEs at impact parameter 50 kpc < R_\_ < 130 kpc 
per quasar is A^lae = 0.34 for our conservative lower 
limit WLya > 50A and TVlae = 0.029 for our formal 
lower limit Wi^ya > 173 A. Note that this calculation is 
conservative and would tend to over-predict the number 
of LAEs near quasars for two reasons. First, we used the 
quasar-LBG correlation function, but LBGs are known to 
cluster slightly stronger than LAEs. Second, quasar clus- 
tering varies rat her strongly with reds hift (Shen e t al.l 
I2OO8D and so the lTrainor Steidell ([2012 ) quasar-galaxy 
clustering measurements at z ^ 2.7 probably overesti- 
mate the clustering around z ^ 2 quasars. 

But given that we expect A^lae sources in the annu- 
lar region between i?_L,min and i?±,max around a given 
quasar, what is the probability that one such source 
should land on our slit given our observed sample? As 
we used our median luminosity limit above (i.e. the 
limit in SDSSJ 0856+1158 scaled up by three) to com- 
pute the number density of LAEs tilae, we assume that 
we could have detected a Li^ya > 6.4 x lO'^^ergs"-^ in 
15/29 sources in our sample. Each longslit observation 
covers 3% of the annular region around a given quasar, 
hence the total number of LAEs expected in our sample 
is 15 X 0.03 X A/'lae, which equals 0.15 LAEs for our con- 
servative EW limit, and 0.013 for our formal equivalent 
width limit. Although our calculation is crude, and does 
not fully take into account the sensitivity variations in 
our sample, it seems rather unlikely that the clustering 
of LAEs around the quasar can trivially explain our dis- 
covery of a relatively high-EW LAE at R± — 131 kpc 
from the f/g quasar in SDSSJ 0856+1158. 

5.3.4. Is SDSSJ 0856+1158 Powered by Fluorescence? 

It is intriguing that the properties of the LAE in 
SDSSJ 0856+1158 resemble those of th e candidate fiu- 
orescing source closest to the quasar in Cantalupo et al.l 
(|2012r ). Namely the two sources are at nearly identi- 
cal impact parameter a nd our LAE is o n ly a a factor 
of four fainter then the ICantalupo et al.l (|2012l ) candi- 
date. Is the luminosity of the SDSSJ 0856+1158 LAE 
consistent with Lya mirror fiuorescence? Taking the dis- 
tance from the quasar to be the impact parameter Rj_, 
we compare the mirror surface brightness (see eqn. [T2|) 
to the observed LLya, to deduce the area, and hence size 
of the emitting 'mirror'. For a spherical source the im- 
plied diameter is c/lae = 2.5 kpc which subtends an an- 
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gle 6>LAE = 0.29". This small size is consistent with our 
assessment that the LAE in SDSSJ 0856+1158 is likely 
spatially unresolved, although we reiterate that the data 
are noisy. Note also that our inferred size depends on 
two unknown quantities, namely the geometric reduction 
factor /gm (we assumed /gm = 0.5) and the line-of-sight 
distance to the quasar R\\ (we assumed R\\ = 0), and 

scales as c^lae oc /gm^^[l + {R\\/ R±)'^Y^'^ . It is thus 
quite plausible that the source could be a factor of ^ 3 
larger and hence marginally resolved. We conclude that, 
within these uncertainties, the luminosity of the LAE in 
SDSSJ 0856+1158 is consistent with the mirror fluores- 
cence interpretation, given our observational constraints 
on its distance, size, and the ionizing luminosity of the 
foreground quasar. 

Again, how likely is it that one such fluorescing source 
should land on our slit given our obser ved sample? Here 
we ca n only speculate, but the single iCantalupo et al.l 
(|2012f ) fluorescing source provides circumstantial evi- 
dence that the number of such sources expected within 
the annular region considered is approximately unity 
^LAE ^ 1- Analogous to our argument at the end of 
the previous section, assuming we had the sensitivity 
to detect such a source in roughly half of our sample, 
and that our slit covers 3% of the annular region, we 
get 15 X 0.03 X TVlae or 0.46. This is a factor of - 3 
larger than the number of LAEs expected from cluster- 
ing around the quasar, and we hence favor the fluores- 
cence interpretation, although we emphasize that our es- 
timates, particularly for the fluorescence, are very crude. 

The two physical scenarios, clustering and fluorescence, 
can be distinguished by obtaining sensitive constraints on 
the rest-frame EW. For instance, if we adopt our formal 
limit on the EW in SDSSJ 0856+1158 W^ya > 173 A, 
then the probability of finding one such LAE in our 
sample due to quasar-LAE clustering is ^ 1%, which 
would make a stronger case for fluorescence. Further- 
more, if it were established that Wi^ya > 240A, then 
this would be very compelling evidence for fluorescence 
because normal stellar populations are not expected to 
produce EWs above this value. Given our current rel- 
atively weak constraints on the EW, it is not possi- 
ble to make a convincing statement about whether the 
LAE in SDSSJ 0856+1158 results from clustering or flu- 
orescence. Deep broad band imaging of the LAE in 
SDSSJ 0856+1158 could be used to obtain a much more 
sensitive constraint on the continuum, distinguishing be- 
tween these two scenarios. 

6. SUMMARY AND CONCLUDING REMARKS 

In this paper we introduced a novel technique to si- 
multaneously analyze the absorption line and emission 
line properties of CGM and IGM gas around a quasar 
using close projected quasar pairs. We analyze the 
absorption spectra of the b/g quasars in a sample of 
68 projected quasar pairs, and identify 29 systems for 
which the spectrum of the b/g quasar shows evidence for 
an LLS coincident with the redshift of the f/g quasar. 
The presence of self-shielding circumgalactic gas, pin- 
pointed by the b/g sightline, makes these 29 objects 
prime candidates for the detection of Lya emission from 
the CGM, either from fluorescent recombinations pow- 
ered by the quasar, pure resonant scattering of Lya pho- 



tons emitted by the quasar, or Lya cooling radiation 
from T ~ 10^ K circumgalactic gas. We perform a sys- 
tematic, slit-spectroscopic survey for extended Lya emis- 
sion in the vicinity {R± < 600 kpc) of these 29 z ~ 2 
quasars, achieving typical la surface-brightness limits of 
SBLya ^3 X 10~^^ ergs~^ cm~^ arcsec"^. The primary 
results of this survey and our analysis are: 

• Fluorescent Lya emission with properties expected 
in the 'mirror' approximation, is not detected at 
the locations of any of the 29 optically thick ab- 
sorbers pinpointed by the b/g quasar sight lines, al- 
though we had the sensitivity to detect this sig- 
nal in 24/29 sources at > lOa significance. We 
thus seriously question the realit y of the purported 
detection of this phenomenon by lAdelberger et al.l 
(|2006h . 

• We measure the covering factor of optically thick 
absorption in the quasar CGM to be /c ^ 0.5 on 
R± < 200 kpc scales from our sample of 29 ab- 
srobers detected in a parent sample of 68 pair sight- 
lines. This high covering factor implies a level of 
large scale ~ 100 kpc diffuse Lya emission that 
would have been easily detected for 28/29 sight- 
lines we searched for emission. Despite this fact, 
we only detected large scale Lya emission in a sin- 
gle system. 

• Both the lack of mirror fluorescence at the discrete 
location of absorption in the b/g sightline, and 
the lack of diffuse fluorescent emission extended 
along the slit at R_\_ ^ 100 kpc can be explained if 
the gas absorbing b/g sightlines is shadowed from 
the f/g quasar ionizing radiation, due to the ob- 
scuration which is expected from unified models 
of AGN. This same obscuration argument (previ- 
ously put forth in QPQl and QPQ2), also explains 
the anisotropic distribution of absorbers around 
f/g quasars (QPQ2). That is, the very high in- 
cidence of optically thick absorption observed in 
b/g sightlines fc ^ 0.5, compared to the much 
lower incidence of strong absorption along the line- 
of-sight. In this context, we crudely estimate the 
average opening angle of our f/g quasar UV emis- 
sion to be 1] < 27r which implies an obscured frac- 
tion /obscured > 0.5, cousistcut with valucs deduced 
from mult i- wavelength studies of AGN. 

• Extended Lya fuzz is detected on small-scales 
R± < 50 kpc at levels of SBLya — 0.5 — 
5 X 10~^^ergs~^ cm~^ arcsec"^ in 10/29 objects, 
or 34% of our sample. Given that our single slit ori- 
entation only covers a fraction of the area around 
the quasar, the true detection rate is expected to 
be significantly higher ~ 50 — 70%, as deduced 
by other recent studies. This Lya fuzz is far too 
faint to be arising from the high covering factor 
fc ^ 0.5 optically thick gas we frequently detect in 
b/g sightlines. If the emission arises from optically 
thick gas, the covering factor has to be very low 
fc ^ 0.01, which could arise from a population of 
rare clouds with very high density sufficient to self- 
shield against the intense quasar radiation; however 
this interpretation is not unique. This emission 
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is likely the high-redshift analog of the extended 
emission-line regions (EELRs) commonly observed 
around low-redshift {z < 0.5) quasars. 

• A Lya-emitter with high-equivalent width Wi^ya > 
50 A (rest-frame) and a luminosity Li^ya = 2.1 ± 
0.32 X lO^^ergs"^ is detected at a distance of 
134 kpc from the f/g quasar in SDSSJ 0856+1158. 
Assuming the Lya is powered by star-formation, 
we estimate the probability of discovering such an 
LAE in our survey, including the enhancement due 
to clustering around the quasar, to be small < 15%. 
This source has properties comparable to the fluo- 
rescent emission candidates recently uncovered by 
iCantal upo et al.l (|2Q^ 12). and their discovery of flu- 
orescent LAE at about the same impact parameter 
from the quasar, provides circumstantial evidence 
that the LAE in SDSSJ 0856+1158 is also powered 
by fluorescence. 

In the next few years, new sensitive optical integral- 
field spectrometers will be coming online on large- 
aperture telescope s, such a s the Multi Unit Spectroscopic 
Explorer (MUSE; |Baconet al. 201 QD and the K eck Cos- 
mic Web Imager fKCWL Martin et al.ll2010f ). These 
image slicing integral-field units (IFUs) will enable 3-d 
spectroscopy over ~ 1' x 1' fields-of-view, comparable 
to the linear angular dimension we mapped out around 
quasars with our l.O^'-wide slit observations. The multi- 
plexing gain in speed of these instruments is thus a factor 
~ 60 for mapping out emission from the quasar CGM. 
It is anticipated that these instruments will achieve SB 
limits SBLyo; — 3 x 10~^^ ergs~^ cm~^ arcsec"^ in ~ 50 
hour integrations. 

In light of the the results in this paper, it is inter- 
esting to speculate about what these instruments might 
detect at these unprecedented sensitivity levels. Con- 
sider a quasar with log^g ^i^Ly« = ^1, which is compa- 
rable to the brightest sources among the 29 we con- 
sidered, corresponding to a typical i-band magnitude 
i = 17.7, or about 1.6 magnitudes brighter than our 
median source. Based on the statistics of absorp- 
tion in background sightlines (e.g. see Figure [2] and 
QPQ1,QPQ2,QPQ5) and our detailed analysis of the ab- 
sorption in a single system (QPQ3), our favored model 
for the properties of the quasar CGM is (i? , nn , Nn , /c) 
= (100 kpc,0.1cm-3,102°cm-2,0.5). If illuminated by 
the quasar, such clouds will be highly ionized and opti- 
cally thin to ionizing radiation, but optically thick to 
Lya photons. The predicted level of optically thin flu- 
orescence is SBLya = 3.8 X 10~^^ ergs~^ cm~^ arcsec"^ 
and scales as the product nuNn. Pure Lya scattering 
gives SBLya = 2.7 x 10~^^ ergs~^ cm~^ arcsec"^. The 
Lya scattering would be easily detectable by these future 
instruments for long integrations, and even the fainter 
optically thin emission would also likely be detectable if 
one averages spatially. We thus conclude that our obser- 
vations and calculations imply that future deep integra- 
tions by MUSE and KCWI are expected to uncover glow- 
ing Lya emission on R ^ 100 kpc scales around nearly 
every bright quasar, allowing one to directly image the 
CGM. The level and kinematics of this emission, when 
combined with absorption line information, will provide 
new information about properties of the quasar CGM, 



and a new window in studying the physics of galaxy for- 
mation in the massive halos hosting quasars. 
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Figure 5. Two dimensional spectra plotted as x-maps for our entire projected quasar pair sample. The bottom and middle panels show 
Xsky (sky-subtracted only) and Xsky+PSF (sky and PSF subtracted), respectively, as defined in 21 the absence of extended emission, 
the distribution of pixel values in the the Xsky+PSF (middle panels) images should be Gaussian with unit variance (see Figure [3]). All 
images are displayed with a linear stretch ranging from —6a to 6a as indicated by the central color bar. Upper panels show smoothed maps 
Xsmth (see which are helpful for identifying extended emission. The stacked images were smoothed with a symmetric Gaussian kernel 
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THE CGM in ABSORPTION & EMISSION 



-5000 



Av ^(km/s) 



5000 



30 
20 
10 

o 

Q) 

^ 

CO 

-10 
-201- 



30 
20 
10 

o 

Q) 

g 

CO 

-10 
-20 

30 
20 
10 

o 

CD 

o U 
CO 

-10 
-20 



L 0841+3921' 'i • ^ 

• . -J. ». . 


■ ' ' - 


i- k' - 1-1 • 


^ *^ 

. 

• " • * : 


|--^ ■ 1 ■ ■ ■ ■ H-' ■ 




r^-----^ T-:--t:-| ----- 
1 ■ ■ 1 ■ ■ ■ ■ 1 ■ ■ ■ 


'J 















300 
200 
100 

-100 

-200 
300 

200 

100 

0^ 

-100 

-200 
300 

200 

100 

-100 
-200 



-6000 -4000 -2000 



-5000 



5000 



Av 



-6000 



Av 

-4000 -2000 



(km/s) 
(km/s) 

2000 4000 



20 
10 

o 

CD 

o U 

CO 

-10 

20 
10 

o 

^ 

CO 

-10 

20 
10 







-10 



'0924+3929 " 



-i — . — . — . — I — . — . — . — r 



6000 
200 



6 
5 
4 
3 
2 
1 

-1 

100 -2 
-3 



I 



-4 

,-5 
-6 



0^ 

-100 
200 

100 

0^ 

-100 
200 

100 

0^30 
-100 



I 



-6000 



-4000 -2000 2000 4000 
A V (km/s) 



6000 



20 
10 

o 

^ 
CO 

-10 
-20 
20 
10 

o 

1 

CO 

-10 
-20 
20 
10 

o 

CO 

-10 
-20 




' ■ ■• ■ I ' 



200 
100 

-100 

-200 
200 

100 
-100 



-200 
1200 



I ■ ■ ■ I ■ ■ ■ I 



100 

0^30 
-100 
-200 



-6000 -4000 -2000 2000 4000 
A V (km/s) 



6000 



A V 

6000 -4000 -2000 
I ' ' 



(km/s) 
2000 4000 
■ ■ I ■ ■ ■ I ■ 



20 i- 
10 1- 



co 

-10^ 
-20 i- 



20 
10 

o 

^ 0^ 

CO 

-10 
-20 



CD 



20 ■- 

loE- 



CO 

-10^ 
-20 i- 



0938+5317 " 



6000 



200 

100 

0^ 

-100 

-200 
200 

100 

0^ 

-100 

-200 
200 

100 

0^ 

-100 

-200 



-6000 -4000 -2000 2000 4000 
A V (km/s) 



6000 



Figure 5. Continued. 



30 

-6000 -4000 



HENNAWI & PROCHASKA 



A V (km/s) 
-2000 2000 4000 



20 
10 

U 



-10 
20^ 
10 

D 

3 

-10 
20 
10 



-10 



1006+4804 



0) 

55 [ i w nitm i im ainmm.m i * 



6000 
"^200 
-i 150 

100 
-Iso 

-i -50 

-: 200 
^ 150 

100 

-i 50 
-i -50 

-i -100 

-: 200 
^ 150 
100 

;50 

; -50 
; -100 



-6000 -4000 



-2000 2000 4000 
A V (km/s) 



A V (km/s) 
-10000 -5000 5000 



10000 



40 



-40 



-40 



. 1 1 1 1 1 1 1 1 1 
: 1041+5630 


. . 1 . . . 1 1 . 


. 1 . 1 L ■ 




* * 




_ " f 






^ 1 1 1 1 • 1 1 1 1 1 








1 1 1 1 1 1 1 1 1 


1 1 1 1 1 ■ 














'. 1 .... 1 .... 1 .... 1 .... 1 . ' 














'. 1 .... 1 .... 1 ... . 1 ... . 1 . ' 



6000 

400 

200 

^ 

-200 

-400 

-600 



6 
5 
4 
3 
2 
1 

-1 



-4 



I 



400 

-5 

200 _6 

^ 

-200 

-400 

-600 

400 

200 

? 

-200 

-400 

-600 



I 



A V (km/s) 

-6000 -4000 -2000 2000 4000 6000 



20 
10 

o 

Q) 

o ^ 

CO 

-10 
-20 
20 
10 

o 

Q) 

o ^ 

CO 

-10 
-20 
20 
10 


CO 

-10 
-20 



CD 



i 1025+5820 


■ 1 ■ 1 ■ ■ 


' 1 J 








M ■ ■ ■ 1 ■ ■ ■ 1 ■ ■ ■ 1 ■ ■ ■ 








1 ■ ■ ■ 1 ■ ■ 












h ■ ■ ■ 1 ■ ■ ■ 1 ■ ■ ■ 1 ■ ■ ■ 








1 ■ ■ ■ 1 ■ ■ 












■ ■ 1 . ; . .1 . ■ . . 1 . ; .1 . . . 


. . . 1 . . .. 


.■ 1 ■ 



200 
100 

-100 

-200 
200 

100 
-100 



-200 
200 



100 

-100 
-200 



-6000 -4000 -2000 2000 4000 6000 
A V (km/s) 
A V (km/s) 



-6000 -4000 -2000 



2000 4000 6000 



20 
10 

o 

CD 

CO n 

o u 

CO 

-10 

-20 

-30 

-40 
20 

10 

o 

CD 

CO n 

CO 

-10 

-20 

-30 

-40 
20 

10 

o 

CD 

CO n 

CO 

-10 
-20 
-30 
-40 



■ I ■ ■ ■ I ■ ■ ■ I ■ ■ ■ I ■ ■ I ■ ■ ■ I ■ ■ ■ I 
1045+435t,^» I ^^.^ # • % - - 



I ■ ■ ■ I ■ ■ ■ I 



I ■ ■ ■ I ■ ■ ■ I 



100 

^? 
-100 

-200 

-300 

100 

^? 
-100 

-200 

-300 

100 

^? 
-100 

-200 

-300 



-10000 -5000 5000 10000 

A V (km/s) 



-6000 -4000 -2000 2000 4000 6000 
A V (km/s) 



Figure 5. Continued. 



THE CGM in ABSORPTION & EMISSION 



Av 

-6000 -4000 -2000 



lOO 4000 6000 



-4000 



A V (km/s) 
-2000 2000 



20 
10 

o 

Q) 

^ 

CO 

-10 
-20 

20 
10 

o 

Q) 

^ 

CO 

-10 
-20 

20 
10 

o 

CD 

g 

CO 

-10 
-20 



1145+0322 . 



: 200 
100 

-100 
-200 
^200 
100 

-100 
-200 
3200 
100 

io^ 

-100 
^-200 



-6000 -4000 -2000 2000 4000 6000 
A V (km/s) 
A V (km/s) 

-6000 -4000 -2000 2000 4000 6000 



20 
10 

o 

CD 

o u 
CO 

-10 
-20 
20 
10 

o 

CD 

CO n 

CO 

-10 
-20 
20 
10 

o 

CD 

g 

CO 

-10 
-20 



1306+6158 " " ' 



-1- 



4 •* 



6 

5 
4 
3 
2 
1 


200 _^ 

-2 

-3 

-4 

-5 
-6 



I 



100 
-100 



200 
100 
0^ 
-100 

200 
100 
0^ 
-100 



-6000 -4000 -2000 2000 
A V (km/s) 



Figure 5. Continued. 



4000 6000 



20 
10 

10 
20 

20 
10 


-10 
-20 
20 
10 

o 

Q) 

^ 
10 
20 



CO 



CO 



1204+0221" 



I ■ ■ ■ I 



MfMM 



-4000 -2000 2000 

A V (km/s) 



4000 
3 200 

100 

0^ 

-100 

-200 
200 

100 

0^ 

-100 

-200 
200 

100 

0^ 

-100 

-200 
4000 



A V 

-6000 -4000 -2000 
20 i 



000 4000 



10 

o 

CD 

o u 

CO 

-10 
-20 
20 
10 

o 

CD 

o u 

CO 

-10 
-20 
20 
10 




-10 
-20 



CD 



1358+2738 " 



-r 



-r 



I ■ ■ ' . I ■ ■ i. i' . ' ■ I ■ ■ ■ ■ I 



6000 
"n200 



100 

0^ 

-100 

-200 
^200 

100 

0^ 

-100 



-200 
200 



100 



mmi p 



-6000 -4000 -2000 2000 4000 
A V (km/s) 



-100 



-200 
6000 



32 HENNAWI & PROCHASKA 



Av (km/s) Av (km/s) 

-6000 -4000 -2000 2000 4000 6000 -4000 -2000 2000 4000 




Av (km/s) _4000 -2000 2000 4000 

A V (km/s) 



Figure 5. Continued. 



Av 

-6000 -4000 -2000 



20 r 
10 I- 

o 

^ oh 

CO 

-10 i- ' 



-20 r 



20 : 

10 \ 



3 oH* 

B 

-10 
-20 



1442+0137 



THE CGM in ABSORPTION & EMISSION 

(km/s) 

2000 4000 6000 

200 



Av 

-6000 -4000 -2000 
20 i 



(km/s^^ 



^100 

^-100 

-i -200 
200 
100 

-100 

-200 

200 

100 

0^ 

-100 

-200 



0) 



-6000 -4000 -2000 2000 4000 6000 
A V (km/s) 
A V (km/s) 



-6000 -4000 -2000 



2000 4000 6000 



10 



Y> 



-10 



-20 



10 



-10 



-20 



10 



1508+3635 



100 
50 

-AO J3 -; 



WO 



-50 

-100 

-150 

-200 
100 

50 

0^ 

-50 

-100 

-150 

-200 
100 

50 



6 
5 
4 

3|- 
2 

1 h 


-1 h 
2 

3L 



10 

o 

0) 

^ 
CO 

-10 
-20 
20 
10 


CO 

-10 
-20 
20 
10 

o 

Q) 

^ 

CO 

-10 



1444+3113 



:000 4000 6000 
200 



100 
0^ 
-100 

200 
100 
0^ 
-100 

200 
100 
0^ 
-100 



-20 

-6000 -4000 -2000 2000 4000 6000 
A V (km/s) 



Av 



-6000 -4000 -2000 
-I- 



(km/s) 

2000 4000 6000 
— 1 — I — I — I — I — I — |— 



I 



CO 



-50 
-100 
-150 
-200 



-6000 -4000 



-2000 2000 4000 6000 
A V (km/s) 



20^ 
10 


-lOt- 
-20 f- 
20 

lot-. 

o 

CD 

^ oE- 

CO 

-10 
-20 
20^ 
10 


CO 

-20 



CD 



1553+1921 



- 200 

-i 100 

- -100 

^ -200 
200 

-i 100 
. ^-100 



-200 
200 



- 100 

■ ^ 

-. -100 
-i -200 



-6000 -4000 -2000 2000 4000 6000 
A V (km/s) 



Figure 5. Continued. 



34 



HENNAWI & PROCHASKA 



A V (km/s) 
-4000 -2000 2000 4000 
-r 




I ■ J ■ I 



200 



100 



-100 



-200 
200 



100 



op 



-100 



-200 
200 



100 



i 



-100 



-200 



-4000 -2000 2000 4000 
A V (km/s) 



Figure 5. Continued. 
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APPENDIX 
THE ATTENUATION OF Lya BY DUST 

In most, if not all, astrophysical environments that contain cool gas and metals, a fraction of the refractory elements 
(e.g. Fe, Ni, Si) have condensed to form dust grains. For the treatment of Lya emission considered in this paper, dust 
is an important factor. Because we detect absorption from gas which is optically thick in the Lya line, Lya photons 
will be resonantly trapped, and can only escape the system through diffusion, both in frequency and spatially. If the 
resonant photons must traverse a much larger pathlength to escape, this could result in a significant increase in the 
extinction of resonant lines compared to the continuum extinction at the same wavelength. Because the focus of this 
paper is Lya emission from the cold gas detected in absorption, we calculate the magnitude of the attenuation of Lya 
by dust grains given our expectations for the physical conditions of the gas in absorbing clouds. Our estimates closely 
follow the discussion in §15.7 of lDraine (2011). 

Consider an ionized spherical cloud with total hydrogen column density A^h, neutral column density A^hi, neutral 
fraction xhi = A^hi/A^h, volume density nn, and Doppler broadening parameter b. We define the radius of the cloud to 
be Vc = 3A^H/4nH, where we have used the fact that the average absorption pathlength (/c) = A^h/^h through a sphere 
of radius Vc equals (/c) = 4/3rc. In what follows we assume, that the cloud has a uniform ionization structure, whereas 
in reality such clouds are optically thick to ionizing photons and will be self-shielding. Uniform ionization structure 
is likely to be a good approximation in a cooling radiation scenario where Lya photons are generated throughout the 
cloud, as the bulk of the Lya photons will be generated in the regions of the cloud which have properties close to the 
mean. For photoionization of self-shielding clouds, this approximation is likely to over-estimate the dust attenuation, 
because most recombinations are generated in a thin self-shielding layer at the surface of the cloud where the neutral 
fraction (and hence Lya optical depth) is much lower than the average value in the cloud. 

Following .Draine (2011), we can write the mean square displacement of a Lya photon from the point of emission as 

{r') - IlIn! (A1) 

Lo ^ (A2) 
nmcTw 

(j^ = 6.48 X lO-^^cm-2 ( (A3) 

V20kms"V 

and A^s is the total number of resonant scatterings and is the H I absorption cross-section. These expressions are 
valid in the Lorentzian wings of the Voigt profile where the bulk of the spatial diffusion takes place. A reasonable 
criterion for escape is that (r^) = r^, from which we can determine the total number of scatterings via eqn. ()Aip 



where 



1 /2 

This modest number of scatterings is sufficient to frequency-shift a typical photon to Au/u ^ Ns b/c, such that our 
assumption of being in the damping wings of the Voigt profile is valid. 

To determine the amount of attenuation by d ust, w e must first estimate the total pathlength traveled by an escaping 
photon, which can be written as 5 ~ l/2LoN^ (jDrainell2011i) . or 

JV„, >. 



Thus given the cloud properties we have assumed, the pathlength and hence total optical depth through the cloud for 
resonantly scattered photons is a factor of 2 larger than it is for continuum photons. 

The amount of extinction corresponding to this pathlength depends on the amount of dust associated with the 
gas, which is highly uncertain. In quasar absorption line research, there is strong evidence that refractory elements 
in DLAs are depleted onto dust grains albe it at modest levels compared to the Galactic ISM (jPettini et al.l [1993 : 
iProchaska fc Wolfell200TI: IVladilo et aII[20TTI ). The effects appear to be strongest for the systems with highest metal- 
licity (e.g. lKaplan et al"]l2010f ). For LLSs where the gas is predominantly ionized, there has been essentially no empirical 
assessment of dust. One expects that the physical conditions - higher UV radiation field, higher temperature, lower 
density - l ead to lower dust-to-gas ratios. Nevertheless, dust is apparent i n highly ionized reg ions of the local uni- 
verse (e.g. iDraindHoilf ) and also the diffuse environment of galactic halos (jMenard et al.ll2010f ). The most extreme 
assumption is to assume that all of the gas contains dust, and that the dust-to-gas ratio corresponds to that of the 
Milky Way, which is equivalent to assuming a solar metallicity and a high depletion factor. This is conservative, 
because DLAs are observed to have lower depletions than the Milky Way ISM and it is rather unlikely that all of gas 
in the CGM of quasars has solar enrichment as we assume. For example, our detailed analysis of SDSSJ1204+0221 in 
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(|PrQchaska fc Hennawill2QQ"9l ) showed met allicities in the ran^e .3 — 1.6 Z©, and our expectation is that a more typical 
metaUicity of the quasar CGM ~ 0.1 Zq (jProchaska et al.ll2012L Prochaska et al., in prep.). 

The continuum extinction at A = 1216A can be written ^1216 = 0.18(i^v/3.1)(A^H/10^^ cm~^) m ag, co rresponding 
to a visual extinction of Ay = 0.05, where we take Al2ie/ Ay = 3.3 based on the extinction curves of (jFitzpatrick 1999) 
for Ry = 3.1. In the absence of resonant effects, 12 16 A photons would be attenuated according to the continuum 
extinction ^1216- Combining this extinction with the increased resonant pathelngth in eqn. ()A5p . we can write the 
total resonant extinction that Lya photons suffer as they ecape the cloud as 

Once Lya photons escape an individual cloud, they propagate freely until they encounter other clouds, where they 
might undergo additional resonant scattering. But the probability of this occurring is small, since it requires both 
spatial and spectral overlap. If the spatial covering factor of clouds is /c, and if bulk motions of the clouds in the 
quasar halo follow a gaussian with line-of-sight dispersion a, then the probability of additional resonant scattering is 

^ fcNs^'^b/V2a - 0.07 for typical numbers {N^ = 15, fc = 0.50, a = 400km s"\ and b = 20km s"^), thus additional 
scattering would amount to < 10% corrections to the total extinction, which is already highly uncertain because of 
the unknown dust-to-gas ratio. However, in traversing additional material Lya photons will still encounter continuum 
extinction, so we can heuristically write the total extinction as A1216 ^ ^me + /c^me- 

Thus given the observed properties of the quasar CGM (A^hi ^ 10^^ cm~^; xm ~ 0.1; b ~ 20 km s~ ), and under the 
most conservative assumptions regarding the amount of dust, namely a Milky Way dust-to-gas ratio (corresponding 
to solar metallicity and high depletion), we conclude that the extinction of Lya photons would be in the range of 
^1216 ^ 0.3 — 0.6, resulting in at most a factor of ~ 2 reduction in the total Lya surface brightness. This extinction 
scales linearly with the metallicity, and as mentioned previously, it is unlikely that all gas in the quasar CGM has solar 
enrichment as we have assumed. Furthermore, our estimate of the extinction due to resonance effects is appropriate 
for Lya photons generated throughout the volume of an emitting cloud, as expected in a cooling radiation scenario, 
whereas it is likely overestimates the attenuation of flourescent emission powered by photoionization. 

Finally, it is worth asking whether the continuum extinction from CGM absorbers would be detectable via the 
reddening of the b/g quasar spectra. C onsider for example the absorber in SDSSJ1204+0221 which we analayzed 
in detail in (jProchaska fc Hennawil [20091 ) . Our photoionization modeling gave a higher total gas column density of 
A^H = 10^^"^ cm~^, which implies a continuum extinction of ^^216 ~ 0.57, for solar metallicity and a Milky Way dust to 
gas ratio. This reddening would not be distinguishable from either the SDSS photometry or our optical spectra, since 
the relative reddening ^^216 — ^2500 = O-^^ between rest frame wavelengths A = 1216A and 2500A (corresponding to 
the observed frame 4000 — 8000A covered by the optical spectra and photometry) is small, and cannot be distinguished 
from intrinsic variations in the quasar continuum slopes. 
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